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Abstract. In standard nonrelativistic quantum mechanics the expectation of the 
energy is a conserved quantity. It is possible to extend the dynamical law associated 
with the evolution of a quantum state consistently to include a nonlinear stochastic 
component, while respecting the conservation law. According to the dynamics thus 
obtained, referred to as the energy-based stochastic Schrodinger equation, an arbitrary 
initial state collapses spontaneously to one of the energy eigenstates, thus describing 
the phenomenon of quantum state reduction. In this article, two such models are 
investigated: one that achieves state reduction in infinite time, and the other in 
finite time. The properties of the associated energy expectation process and the 
energy variance process are worked out in detail. By use of a novel application of a 
nonlinear filtering method, closed-form solutions — algebraic in character and involving 
no integration — are obtained for both these models. In each case, the solution is 
expressed in terms of a random variable representing the terminal energy of the system, 
and an independent noise process. With these solutions at hand it is possible to 
simulate explicitly the dynamics of the quantum states of complicated physical systems. 
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1. Introduction 



The idea that the standard Schrodinger equation of nonrelativistic quantum mechanics 
should be extended to take the form of a stochastic differential equation on Hilbert 
space has been investigated extensively as a mathematically viable approach to the 
measurement problem in quantum mechanics. Indeed, there is now a substantial body 
of literature on the theory of spontaneous state vector reduction, and a number of 
different models have been proposed that fall into this category. See, e.g., [SI 13 EH Ell 
for overviews of this area and relevant references. 

This paper is concerned with the so-called energy-based stochastic extension of the 
Schrodinger equation, which has the special status of being the simplest such extension 
that is in principle applicable to any nonrelativistic quantum system. The physical setup 
can be described briefly as follows. We consider an isolated quantum system for which 
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the Hamiltonian H has a discrete spectrum {-£?i}i=i ) 2,...,jv- We assume that initially the 
system is in a pure state represented by the state vector l^o)- The situation where the 
initial state is mixed can also be considered (see, e.g., jlj), but for simplicity we confine 
the discussion to the case of an initially pure state in this paper. For convenience we 
set (00 1 Vto ) = 1. For each value of i we let |0j) denote the Liiders state associated with 
the energy level E^. More specifically, let us write Ilj for the projection operator onto 
the Hilbert subspace of states for which the energy has the value E^. We allow for the 
possibility that the energy may be degenerate. In that case, we write rii for the number 
of linearly independent state vectors with energy Then the Hamiltonian takes the 
form H = J2i Eifli, and we define the Liiders states by setting 

10.) = ft#0) ■ (1) 
(^IWo) 1 / 2 

We note that H\<j>i) = Ei\<f>i) and (0i|0t) = 1. According to the von Neumann-Liiders 
projection postulate j22E0I, if the system is initially in the pure state |-0o) an d if the 
outcome of a measurement of the energy is the eigenvalue Ei, then after the measurement 
the state of the system will be the Liiders state |0j). 

It is an implicit feature of the projection postulate that quantum evolution 
progresses in accordance with the unitary dynamics of the Schrodinger equation up 
to the moment when a measurement is made, at which point the system jumps to a new 
state. The von Neumann-Liiders rule asserting how the jump proceeds is essentially 
ad hoc in nature, despite being plausible from a physical point of view insofar as 
the predicted outcome is concerned. Thus although the projection postulate, in one 
form or another, remains an accepted part of the everyday use of quantum theory in 
practical applications |19j . one has to agree that such a 'cookbook' approach to the 
measurement problem is ultimately unsatisfactory; and this is why over the last five 
decades many attempts have been made to modify the dynamics of standard quantum 
mechanics in such a way that the 'collapse' process can be understood as governed by 
an evolutionary law that operates on a universal basis, rather as does the Schrodinger 
equation in ordinary quantum mechanics. 

In order to ensure consistency with established facts, such a universal evolutionary 
law needs to have the property that for some systems it proceeds in a way that for 
all practical purposes reproduces the dynamics of the Schrodinger equation, whereas for 
other systems the evolution progresses continuously to a terminal state that is consistent 
with the action of the projection postulate. This 'viability' property is satisfied in 
particular by the standard energy-based stochastic Schrodinger equation. In this model 
the Schrodinger equation, which when written in differential form is given by 

d|^> = -iH\ip t )dt, (2) 

is generalised and elevated to the status of a nonlinear stochastic differential equation 
on Hilbert space: 

d\^ t ) = -iH\^)dt - \a 2 {H - H t ) 2 \^ t )dt + \a{H - H t )\^ t )dW t . (3) 
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Here {\i>t)}o<t<oo is the state- vector process, {W t }o<t<oo is the Wiener process, and 
{H t } < t<oc is the energy expectation process, defined by 

Ht - ~mr- ( ' 

The coupling constant a appearing in (jSJ), which has the units a ~ [energy] -1 [time]" 1 / 2 , 
determines the characteristic timescale tr associated with the rate of collapse of the 
wave function induced by (jBJ). This timescale is given typically by an expression of the 
form tr = l/cr 2 Vo, where Vq is the initial value of the process {Vt}o<t<oo for the energy 
variance, which at time t is given by the following expression: 

n = TTTT~\ • 

mm) 

One of the attractive features of the stochastic differential equation © is that 
it provides a more or less completely tractable model for state vector reduction in 
nonrelativistic quantum mechanics. Nevertheless, despite the fact that the mathematical 
properties and physical consequences of (|HJ) and various related processes have been 
studied extensively in the literature [H 13 IU 13 HOI US US HZl HH1 , it has only been 
recently that a general solution of Q has been obtained in terms of an appropriate set 
of freely specifiable random data [TTJ H2] . The aim of this paper is to present a complete 
treatment of the method of solution of the dynamical equation (jSJ). The results are of 
interest both for the new range of numerical and computational techniques they open 
up, as well as for the new methods for model building they provide. 

The article is organised as follows. In §2] we review some basic notions of stochastic 
analysis, including the concepts of nitrations, conditional expectations, martingales, 
supermartingales, and potentials. The material introduced in this section will be 
used throughout the paper. In SI SI and SI we establish various properties of the 
energy expectation process (j3J) and the energy variance process (0), showing that the 
variance process has the 'potential' property; that is to say, its expectation goes to zero 
asymptotically. This allows us to give a precise sense to the notion of state reduction. 
In Sl we then determine the circumstances under which the state vector reduces to one 
of the eigenstates of an observable that is compatible with the Hamiltonian. 

In Sl we address the problem of the origin of the dynamical equation ©. Starting 
from a general stochastic equation for a state vector driven by a single Brownian motion, 
we determine what additional physical assumptions and other simplifying features are 
required in order to obtain (jHJ). We also show, under a suitable 'universality' assumption 
regarding energy conservation, that reduction to lower energy uncertainty is a generic 
feature of the stochastic Schrodinger equation. 

The projection operators for the energy eigenstates constitute a special set of 
observables that commute with the Hamiltonian. The expectation value, with respect 
to the state \ipt), of such a projection operator determines the random process for 
the associated conditional transition probability to that eigenstate. The properties 
of this conditional probability process are studied in SI In SI we digress briefly 
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to investigate the dynamics of the Shannon entropy associated with the system of 
transition probabilities, and show that the Shannon entropy has the property that its 
expectation goes to zero asymptotically. This is contrasted with the behaviour of the 
von Neumann entropy. We also derive an inequality that relates the entropic measure 
and the variance-based measure of energy dispersion. In §101 we study a certain linear 
stochastic differential equation for the state vector, which we call the ancillary equation, 
and verify that a state vector satisfying the ancillary equation, once suitably normalised, 
can be used as a step to obtaining the solution of the nonlinear equation (JHJ). We also 
clarify the relation of our results to earlier work on solutions of (J3J) and related dynamical 
equations, explaining why previously established integral representations for the state 
vector satisfying the stochastic equation should not generally be regarded as explicit 
solutions in the sense that we use the term here. 

Then in §111 and §121 we derive a bona fide explicit solution to (JHJ), making use of a 
nonlinear filtering method. The solution thus obtained is expressed in terms of a simple 
algebraic function of a standard Brownian motion and an independent random variable 
representing the terminal value of the energy. By use of this result, it is possible to 
simulate solutions that represent the evolution of rather complicated quantum systems. 
In Hl'd\ we introduce a technique that allows one to verify that the solution obtained 
in £1111 does indeed give rise to the reduction of the state vector. In ^141 we investigate 
properties of the asymptotic random variable corresponding to the terminal value of the 
energy. While in the nonlinear filtering method used to solve (J3J) we introduce a noise 
term, in §151 we derive the external noise term from the underlying processes specified 
in (J3J). This result justifies the use of the filtering methodology we have employed here. 

We then turn to solve the problem of constructing a collapse model that achieves 
state reduction in finite time. That is, although the standard energy-based collapse 
model achieves a strict collapse in infinite time, with a minimal modification of 
the dynamical equation (JHJ) it is possible to formulate a finite-time collapse model. In 
^161 we introduce such a model. Using the methodology of ^111 we derive an analytical 
expression that we conjecture to give the energy expectation process. The validity of 
this conjecture is established in < J17I In H1S\ we derive the external noise term arising 
in the finite-time collapse model that is used in flYI to obtain the solution. In fT9l 
we demonstrate the fact that the standard energy-based model (JHJ) and the finite-time 
collapse model introduced in ()198|) are related by a nonlinear time-change. That is, 
if we take the model (JHJ) and replace the time variable t by a 'clock' variable defined 
by r(t) = tTjiT — t), where T is a finite positive constant, then in a physical world 
measured by the variable t, the collapse for the new system takes place in finite time 
interval T, since the clock variable r(t) runs from to oo as t runs from to T. 

As a closing remark, in the role of the asymptotic value of the energy, which 
has the interpretation of a hidden variable in the stochastic framework, is discussed. We 
also speculate on whether the energy-based reduction models analysed here suffice as 
such to form a basis for the general description of random phenomena in nonrelativistic 
quantum mechanics. 
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2. Stochastic essentials 

We begin with an overview of the probabilistic framework implicit in the specification 
of the energy-based stochastic Schrodinger equation. The concepts introduced in this 
section are standard in the literature of stochastic analysis, as is also the notation (see, 
e.g., H3 ED EHl HE E3 EH!)- The dynamics of the state vector \ip t ) are defined on a 
probability space (fi, J 7 , P) with filtration {J r t}o<t<oo, with respect to which the process 
{W t }o<t<oo is a standard Brownian motion. Here f2 is the sample space on which T is a 
u-algebra of open sets upon which the probability measure P is defined. Each element 
u) G £1 represents a 'possible outcome of chance'. Each element A G T is an 'event'. 
The measure P assigns a probability ¥(A) to each event A. 

Now we give the relevant definitions in more detail, since these are of interest. Let 
Q be a set, and let JF be a collection of subsets of Q. For any subset A C Q we let 
A c = {uj G Q\u A} denote its complement. Then T is called an algebra of subsets of 
tt if (a) tt G J 7 , (b) A G T implies that A c G J 7 , and (c) A, B E J 7 implies A U 5 G J. 
It follows from these axioms that G J 7 , and that A, B G J 7 implies A R -B G JF. 

The algebraic operations on the elements of J 7 are as follows. The product of two 
elements A, B G T is defined by A - 1? = Af]B, and the sum of two elements A, B G T is 
defined by A + £> = (AU5)fl(An -B) c . It follows that the product and sum operations 
are symmetric and associative, and that (A + B) -C = A- C + B ■ C for any A,B,C G T . 
The underlying field of the algebra T is the minimal subalgebra {f2,0}, which when 
endowed with the same product and sum operations as those defined above satisfies the 
rules of binary arithmetic: 



If J 7 is an algebra of subsets of VL then we say that T is a a-algebra if it has the 
property that whenever {A n } neN G T then U n v4 n G T . That is to say, the union of any 
countable sequence of elements of T is also an element of T . It follows that whenever 
{A n } nm G T we have n n A n G since n n A n = (U n A c n ) c . 

We comment briefly on the distinction between a a-algebra and a topology, since 
the latter is more familiar to physicists than the former. In a topology we axiomatise the 
notion of an open set, and require that the union of any collection of open sets is open, 
and that the intersection of any finite collection of open sets is open. In a a-algebra we 
axiomatise the notion of a measurable set, and require that the union of any countable 
sequence of measurable sets is measurable, and that the intersection of any countable 
sequence of measurable sets is measurable. 

If T is a u-algebra of subsets of a set Q, then we call the pair (Q, T) a measurable 
space. If {Vt,T) is a measurable space, then a probability measure on is a map 

P : T -> [0, 1] satisfying: (a) P(0) = 0, (b) P(O) = 1, and (c) if {A n } neN is a countable 
sequence of disjoint elements of T with union A = U n A n , then F(A) = ^] n P(Ai)- A 
triple (Q, J 7 , P) is called a probability space. 




n ■ n = n, 
n + n = 0. 



(6) 



Quantum noise and stochastic reduction 



6 



The introduction of the concept of a filtration on a probability space allows 
one to formalise the notion that the consequences of the outcome of chance are not 
necessarily revealed at once, but rather may emerge sequentially as time progresses. 
More specifically, a filtration of T is a collection {Tt\ of cr-subalgebras of JF such that 
T s C Tt for all s and t such that < s < t < oo. 

If an event A G T is such that A ^ T t for some given value of t, then we interpret 
this to mean that at time t one can say whether uj G A or not. To put this another 
way, in a filtered probability space each u G f2 corresponds to a possible 'future history'. 
Each element A 6 J 4 then represents a simple yes/no question, the answer to which, 
for any particular future history, will be known for certain by time t. For that reason, 
the nesting T s C T t for s < t gives rise to a notion of causality. 

A real-valued function X : Q i— > R is said to be measurable with respect to the 
cr-algebra .F if for each number x G R the set {w G O : lx(w) = 1} is an element of 
T . Here l x (w) is the indicator function on Q for the set consisting of all uj such that 
X{oj) < x. Thus l x : O h-> {0,1}, and l x (u;) = l{x(w)<x}- If is .F-measurable in 
the sense just discussed, we say that X is a real- valued random variable on (f2,.F, P). 
The probability distribution function F x (x) = P(X < x) is then defined by use of the 
Lebesgue integral: 



More generally, we also consider maps of the form X : i— > F where (Jl, J 7 , P) is a 
probability space and (T, £/) is a measurable space. For example, T could be R n , C n , a 
Hilbert space, or a manifold. In that case we say the random variable X takes values 
in r, and Q can be typically taken to be the so-called Borel cr-algebra generated by the 
open sets of V. Then for any element G G Q we define 



A parametric family {X t }o<t<oo of random variables on (Q, JF, P) is called a random 
process. If a random process {X t } on a probability space (Q, J 7 , P) with filtration {JF t } 
has the property that for each value of t the random variable X t is .F r measurable, then 
we say that {X t } is adapted to the filtration {.Ft}. 

If X is a nonnegative real random variable, its P-expectation (i.e. its expectation 
with respect to the measure P) is defined by the integral 



which may take the value +oo. More generally, when X is not necessarily nonnegative, 
the expectation is defined only when one of the expressions E[X + ] or E[X~] is finite, 
where X + = max(X, 0) and X~ = - min(X, 0), in which case E[X] = E[X + ] - E[X~\. 
A random variable such that E[|A|] = E[X + ] + E[X~] is finite is said to be integrable. 

Now we turn to the definition of conditional expectation. Given a random variable 
X on (fi,jF, P) for which E[X] exists, the conditional expectation E[JT]^4] of X with 




(7) 




(8) 




(9) 
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respect to the cr-subalgebra A C T is defined to be any ^4-measurable random variable 
Y for which E[Y] is defined, such that for any element A e A we have 



If such a random variable exists, then it is unique up to equivalence modulo differences 
on sets of P-measure zero. Thus even if E[X|*4] is not quite unique we refer to it as 
the conditional expectation of X with respect to A. This definition, which at first 
glance appears rather formal and indirect, is nevertheless one of the cornerstones of 
modern probability theory, and is indispensable. We remark that a sufficient condition 
for E[X|.4] to exist is that X should be integrable. 

The following properties of the conditional expectation are often useful in 
calculations: (i) the law of total probability E[E[X|*4J] = E[X]; and (ii) the tower 
property, which says that if A C B C T then E[E[X|£] |„4] = E[X|.4]. The law of total 
probability is a special case of the tower property. 

The conditional expectation operation allows us to introduce the concept of a 
martingale, the stochastic analogue of a conserved quantity. For this purpose we need 
the operation of conditioning with respect to a cr-subalgebra Tt belonging to a filtration 
{J r t}o<t<oo- Intuitively, conditioning with respect to T t means conditioning with respect 
to the information that will become available up to time t. For convenience, we often 
use the abbreviation E t [X] = E[X|jF 4 ] when the choice of filtration can be taken as 
understood. There are situations, however, where more than one filtration may arise 
in the context of a given problem, in which case the more explicit notation is useful. 
The conditional expectation E t [X] satisfies E[E t [X]] = E[X] and E S [E 4 [X]] = E S [X] for 
s <t. We note that if X is ^-measurable, then E([X] = X. 

A real- valued process {X t } is said to be an {Tt}- martingale on the probability space 
(fl, if E[\X t \] < oo for all < t < oo, and E s [X t ] = X s for all < s < t < oo. In 
other words, {X t } is an {jF f }-martingale if it is integrable and if for t > s the conditional 
expectation of X t , given T s , is the value X s of the process at time s. A process {X t } is 
an {T}- supermartingale on (fl,T, P) if E[|X t |] < oo for all t > 0, and E s [X t ] < X s for 
all < s < t < oo. Intuitively, a supermartingale is a process that tends on average, 
at any time, to be nonincreasing. A martingale is a fortiori a supermartingale. The 
martingale convergence theorem (see, e.g., j2H|, theorem 10) states that if {X t } is a 
supermartingale that satisfies 



then there exists a random variable Y such that lim^oo X t = Y almost surely (i.e. 
with probability one), and that E[|y|] < oo. It follows that a positive supermartingale 
necessarily converges to a limit as t goes to infinity. A positive supermartingale {X t } 
with the property that linu^oo E[X] = is called a potential [25] . 

We now review some basic formulae arising in the theory of Ito processes. Let 
(ft, J 7 , P) be a probability space with filtration {^}o<f<oo, an d let {W t }o< t <oo be a 
standard Wiener process adapted to {Tt}- Here by a standard Wiener process (or 




(10) 



sup E[\X t \] < oo 



(11) 



0<t<oo 
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Brownian motion) we mean a continuous process {W t } with the properties that (i) 
{Wt} has independent increments, and that (ii) W t — W s for 0<s<t<ooisa 
Gaussian random variable with mean zero and variance t — s. 

Let {a t }o<t<oo and {&t}o<t<oo be {JF t }-adapted processes such that for any t e [0, oo) 
we have 

/ \a s \ds + [ b 2 s ds < oo (12) 
Jo Jo 

almost surely. Then letting Xo an jF -measurable initial condition, the random variable 
X t defined by the stochastic integral 

X t = X + I a s ds + [ b s dW s (13) 
Jo Jo 

is well defined and ^-measurable for all t, and we call {X t } <t<oo an Ito process (see, 
e.g., [2BJ 1231 f° r the general definition of stochastic integration). In this case we say 
that {Xt} is a real- valued Ito process driven by the one-dimensional Wiener process 
{Wj}. It is straightforward to generalise (|13j) to cases for which both {X t } and {Wt} 
are multidimensional. 

One useful tool of which we make repeated use is Ito's lemma. Suppose {X t } 
is given by (fT3Jl and consider the process {ft}o<t<oo defined by f t = f(X t ,t) where 
/ G C 2,1 (R x ]R + ). Let prime and dot denote differentiation with respect to the first 
and second arguments of f(x,t), respectively. Then Ito's lemma states that 

f(X t ,t)=f{X o ,0)+ [ fXX s ,s)ds+ [ a s f'(X s ,s)ds 

Jo Jo 

+ | f b 2 J"(X s ,s)ds+ [ b s f(X s ,s)dW s . (14) 
Jo Jo 

It is often convenient to express (jlHjl and ()14j) in differential form: thus we write 

dX t = a t dt + b t dW t (15) 

for the 'dynamics' of {X t }, and 

df(X t , t) = (/(X t , t) + a t f'(X t , t) + \b\f\X u t)) dt + b t f'(X t , t)dW t (16) 

for the dynamics of {ft} implied by Ito's lemma. As in ordinary calculus, the differential 
equations of stochastic calculus are essentially formal in character, and always derive 
their meaning from associated integral equations. Thus (fTo^l and (|T6*|) refer back to 
(113)) and (|14j) . Nevertheless, as in ordinary calculus, the manipulation of infinitesimal 
quantities in stochastic calculus can be very powerful as a mathematical technique, and 
can be intuitively very suggestive as well. For example, the so-called Ito product rule 

d(X t Y t ) = Y t dX t + X t dY t + dX t dY t (17) 

is short-hand for the fact that if {X t } is given by (|T3*)l and {Y t } is given analogously, 
but with {a t } and {b t } replaced by {p t } and {qt}, then 

X t Y t = X Y + [ {Y s a s + X sPs + b s q s )ds+ [ (Y s b s + X s q s )dW s . (18) 
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Now consider an Ito process {M t } < t<oo of the form 
M t = M + / b s dW s . 



(19) 



o 



Then a sufficient condition for {M t } to be a martingale is that M should be integrable, 
and that 

' < oo (20) 



E 



bids 







for all t G [0, oo). In that case {M t } is called a square-integrable martingale, and we 
have the identity 

,-t 



E [{M t - M ) 2 ] = E 



bids 



(21) 



More generally, we also have the following relation, valid for t > s > 0, which we call 

the conditional Wiener-Ito isometry: 

(22) 



E s [(M t - M s ) 2 ] = E s I b 2 u du 

J s 



In certain situation we are presented with an equation of the form (|13p. and we 
are told the distribution of X and that the processes {a t } and {b t } are of the form 
at = a(X t , t) and bt = b(X t , t), where a(x, t) and b(x, t) are prescribed functions. In that 
case we have a stochastic differential equation of the form 

dX t = a(X t , t)dt + b(X t , t)dW u (23) 

with initial condition Xq. By a 'solution' of the stochastic differential equation (|2*H|) we 
mean the specification of the probability space (fi,^ 7 , P) with filtration {J-"t}, together 
with an {jF 4 }-adapted Brownian motion and an {jF t }-adapted Ito process {X t } satisfying 
(|23j) along with the given initial condition. 

The extension of these definitions to situations where {X t } and {W t } are multi- 
dimensional is straightforward. It is also appropriate in some circumstances to consider 
processes defined over a finite time horizon t G [0, T], T < oo, for which straightforward 
modifications of the relevant definitions can also be formulated. 



3. Dynamics of the energy process 

Now we are in a position to analyse the dynamics of the energy-based stochastic 
Schrodinger equation (jSJ) in more detail. We shall make the following assumptions 
concerning the dynamics of the state vector: 

(a) the state- vector process {|V't)}o<t<oo takes values in a finite-dimensional complex 
Hilbert space, and is defined on a probability space (fl, J 7 , P) with filtration 
{•T 7 t}o<t< 

(b) {iV^t}} is adapted to {Ft}] and 

(c) (l^t)} satisfies the stochastic differential equation © with the given initial 
condition \ipo). 
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Under these assumptions it is a straightforward exercise in Ito calculus (see, for 
example, to show that (iptl^t) — 1 for all t 6 [0, oo). One is then led to the following 
basic result. 

Proposition 1. The Hamiltonian process {H t } is an {Ft} -martingale, and the variance 
process {V t } is an {J~t} -supermartingale. 

Proof. We need to show that {H t } satisfies 

E s [H t ] = H s , (24) 

and that {V t } satisfies 

KWt] < K, (25) 

where E t [— ] denotes conditional expectation with respect to the cx-algebra Tt- The 
validity of these properties can be established as follows. By an application of Ito's 
lemma to (@J) and (0), we infer that 

dH t = aV t dW t , (26) 

and that 

dV t = -a 2 V t 2 dt + (JK t dW t . (27) 
The process {n t } defined here by 



(28) 



measures the skewness of the energy distribution. Integrating (|2T)jl and (|2*7j) we deduce 
that 

H t = H + a [ V u dW n , (29) 



and that 



Then on account of the relation 



V t = V -a 2 V 2 du + a / K u dW u . 



E, 



b u dW u 



b u dW u 



(30) 



(31) 



o 



that holds for the stochastic integral of any {J-" t }-adapted process {bt} satisfying 

E f b 2 u du < oo, (32) 
Jo 

we deduce the martingale condition (|2^j) from (|2D|) . This follows from the fact that {Vt} 
is bounded. Similarly, it follows as a consequence of (J3T?|) . and the fact that {n t } is 
bounded, that 



E s [V t ] = V s - cr 2 E s 



V 2 du 



which then implies the supermartingale condition (|25|). 



(33) 

<> 
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4. Convergence of the energy variance 

In the case of the Schrodinger equation with a time-independent Hamiltonian, the energy 
process defined by (jH) is constant. This is usually interpreted as the quantum mechanical 
expression of an energy conservation principle. The martingale relation (}24"|) arising 
in the case of the energy-based stochastic Schrodinger equation can be viewed as a 
refinement of this principle. 

The supermartingale property (|25jl satisfied by the variance process is the essence 
of what is meant by a reduction process. In the case of the Schrodinger equation with a 
time-independent Hamiltonian, the variance of the energy is a constant of the motion. 
In other words, not only is the expectation value of the energy fixed, so is the spread. 
On the other hand, the spread of the energy is reduced in the case of the stochastic 
dynamics of (jBJ). In fact, the following result follows as a consequence of equation (fHUJ) . 

Proposition 2. The asymptotic value of {V t } is given by: 

lim E [V t ] = 0. (34) 

Therefore, the variance process is a potential. 

Proof. First we note that if X and Y are integrable random variables, and if X < Y 
almost surely, then E[X] < E[Y]. It follows thus from the supermartingale condition 
PHjl by use of the tower property that if t > u then E[Vt] < E[V U ]. We note that 

r-t 



E 



b u dW u 



o 



(35) 



for any {jF t }-adapted process {b t } satisfying (j3^j) . Since the energy skewness process 
{/tj} is bounded, it therefore follows from (j3UJ) that 

,-t 



E [V t ] =V - a 2 E 



Uo 

t 



V^du 

= V -a 2 I E [V*] du. (36) 
Jo 

Here we have used Fubini's theorem to interchange the expectation and the integration. 
As a consequence, we have the relation 

E[V t ] <V -a 2 [ (E[V u ]) 2 du, (37) 

since E[VJ 2 ] > (E[VJ]) 2 , which follows from Jensen's inequality. Now writing 

v = lim E[V t ], (38) 

let us suppose that v ^ 0. Because E[T4] is a nonnegative, nonincreasing function of 
time, we have v < E[V t ]. It follows from (jHZj) that E[V t ] < V - a 2 v 2 t, which, if v ^ 0, 
implies that E[V t ] vanishes at t — Vo/a 2 v 2 . However, this is incompatible with the 
assumption that v 0; it follows that v = 0, and thus that {V t } is a potential. 
The same conclusion can be reached by a slightly different line of argument. 
Starting with (jHZj), we use the fact that E[V t ] < E[V U ] for t > u to infer that 
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E[T4] < Vq — <x 2 t(E[V^]) 2 , which implies, on account of the positivity of E t [Vf], that 

m\ < v% ( 39 ) 

and hence the claim of the proposition. 

Since V t is nonnegative, Proposition 2 implies that lim^oo V t = almost surely, i.e. 
that reduction to a state of vanishing energy uncertainty occurs with probability one. 



5. Asymptotic properties of the energy 



From the martingale convergence theorem for square-integrable martingales (see ^2} it 
follows that there exists a random variable defined by 



lim Ht 



t— >oo 



(40) 



which represents the terminal value of the energy once the reduction is complete. Thus 
if we write 



(41) 



H + a / V u dW u , 
Jo 

it follows as a consequence of (jSH) that 

H t = E t [H 00 ). (42) 

Thus has the property that it closes the martingale {H t }. It then follows from 
(|4*2*|) that the random variable H t has the interpretation of being the jF r conditional 
expectation of the terminal value of the energy. In particular, we deduce that H = 
Efifoo], which shows that the expectation value of the Hamiltonian in the initial state 
agrees with the expectation of the terminal value of the energy. This result can be 
viewed as a justification for the conventional interpretation of the expectation value of 
the Hamiltonian. 

A similar result can be established in the case of the variance, which we now proceed 
to derive. In particular, writing in the form 



Ef[Vr] = V t — cr E t 



it follows that 



lim Et \V T \ 



v t 



Vfdu 



a z lim Et 



V 2 du 



(43) 



(44) 



Since the variance of the energy is bounded, we can invoke the conditional form of the 
bounded convergence theorem to interchange the order of the limit and the expectation 
on the left-hand side of this equation. It follows from the fact that lim-r^oo Vr = 
almost surely that 



v ± 



a 2 lim Et 



V'fdu 



(45) 
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Now we interchange the order of the limit and the conditional expectation on the right- 
hand side of this equation by using the conditional form of the monotone convergence 
theorem, and we deduce that 



V, 



a 2 E t 



On the other hand, it follows as a consequence of (jSHj) and (f^T^ that 



a 



(46) 



(47) 



(4f 



Equating the results ()46j) and (}48j) we obtain the fundamental relation: 

Proposition 3. Let {iV't}} satisfy (J3J) and write for the asymptotic value of the 
energy martingale {H t }. Then the squared uncertainty of the energy in the state \ipt) is 
given by 

V t = E t [(H 00 -E t [H 00 }) 2 ]. (49) 



V u dW u , 

and therefore, by use of the conditional Wiener-Ito isometry (|22|) . that 
Et [(H^ - H t ) 2 ] = a 2 E t 







"POD 




= a 2 E t 


it v " du 









This relation shows that the random variable V t has the interpretation of being the 
conditional variance of the terminal value of the energy. In particular, Proposition 3 
demonstrates that the initial squared energy uncertainty Vq agrees with the variance 
of the terminal value of the energy. This fact can be viewed as a justification for the 
conventional interpretation of the energy uncertainty. 



6. Asymptotic properties of observables that are compatible with the 
energy 

We now proceed to derive a rather more general result that includes Proposition 3 as a 
special case. Let us suppose G is any observable that commutes with H, and write G t 
and V t G for the mean and variance of G with respect to the random state \ipt)- Thus 
Gt = (G) t and V t G = ((G — G t ) 2 )t, and by use of Ito's lemma we deduce as a consequence 
of © that 

dG t = v-ytdWt, (50) 

and that 

dV t G = -a 2 7 t 2 dt + aStdWt, (51) 

where 

lt ={{G-Gt){H-H t ))t (52) 

and 

6t=((G-Gt) 2 (H-H t ))t. (53) 
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Now we shall show that {G t } is a martingale, and investigate the nature of the 
conditional variance representation admitted by {V t G }. It follows from (joTIj) and (Jo"Tj) 
that 

G t = G + a [ lu dW u (54) 



and that 



V t G = V G - g 1 I ridu + a I 5 u dW u 



tG 



(55) 



Thus, since {74} and {S t } are bounded, we see that {G t } is an {jF t }-martingale and that 
is an {^-"tj-supermartingale. Therefore, by the martingale convergence theorem 
(see there exist random variables Goo and V G such that 



Goo = G + a / -f u dW u 



and 



Vr 



G 



Vr 



G 



a 1 \ 7 ^dw + a \ S u dW u . 







(56) 



(57) 







Taking the conditional expectation of each side of this equation with respect to jF t , we 
deduce that 



00 J 



V ( ? - a 2 E t 



lu du 



+ a / 5 u dW u . 



(58) 



Solving (jHHI) for V G and substituting the result into (jo^|) we see that 



V, 



G 



00J 



E t [V G ]+a 2 E t 



ludW u 



by use of the conditional Wiener-Ito isometry. Making use of the fact that 
Goo — Gt + o 



(59) 



ludW u , (60) 

which follows from ()54|) and ()56|). we obtain the following fundamental relations 
governing the dynamics of {Gt} and {V t }: 

G t = E t [Goo] (61) 

and 

V t G = Et[V G ]+E t [(Goo-Gt) 2 ]. (62) 

Equation (|HT|l shows that the martingale {Gt} closes, and hence the relation 
Go = E[Goo] allows us to identify the initial expectation value Go with the expectation 
of the result obtained for the random variable Goo. Equation (jH^j) then has a natural 
interpretation as a conditional variance relation. If the Hamiltonian has a nondegenerate 
spectrum, then the terminal state is necessarily both an eigenstate of H and G, and V c 



G 

00 
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vanishes. On the other hand, if H has a degenerate spectrum, then the terminal value 
of H will not necessarily be an eigenvalue of G. In that case, the random variable 
is nonvanishing, and takes the value 

Vg = {k\&\4>i) - {4>i\G\<l)iY (63) 
with probability 

tt^KVqI^)! 2 , (64) 

where \<pi) is the normalised Liiders statef corresponding to the eigenvalue Ei of H, 
given the initial state \i[>o), as defined in ^ We recall that when H has a degenerate 
spectrum, the collapse of the wave function induced by (jHJ) necessarily leads to one of 
the Liiders states, as shown, e.g., in 

It is interesting to observe that, while the variance process {V t } associated with 
the Hamiltonian is a potential, the variance {V^ G } for G, although a super martingale, 
is not necessarily a potential unless H has a nondegenerate spectrum. Physically this 
is because a reduction of the energy induces a complete reduction of a compatible 
observable only if the energy spectrum is nondegenerate. 



7. On the generality of the dynamical equation 

Before embarking on an account of our approach to the solution of the stochastic 
differential equation (JHJ), it will be useful to set this stochastic equation in the context 
of a more general family of possible dynamical laws for the state vector process {(V^)}- 
The idea then will be to see what specific additional physical assumptions are needed 
to imply that the dynamics should take the form (j3J). 

We shall assume as before that {\ipt)}o<t<oo is a continuous stochastic process 
defined on a fixed probability space (fi, J 7 , P) with filtration {Tt} taking values in a 
finite-dimensional complex Hilbert space. For the dynamics of (l^t)} we write 

d|V t ) = Mt)dt + aMdW u (65) 

where {W t } is an {jF t }-Brownian motion. We assume that {|^)} is adapted to {Ft}, 
and that so are the operator-valued processes {fit} and {cr t }. We call {fit} and {a t } the 
operator- valued drift and volatility of {IV't)}- 

| See |21]. We remark, incidentally, that the Liiders state has the following geometrical 

interpretation. In the case of a Hilbert space of dimension n + 1, the corresponding space of pure states 
is the complex projective space CP™. If the Hilbert subspace of state vectors of some given energy 
Ei has dimension k + 1, then the corresponding space of pure states of that energy is a projective 
hyperplane D k of dimension k. The complex conjugate of D k is a hyperplane £)"- fe -i f dimension 
n — k — 1. Clearly D k and ]j n - k - 1 & Q n ot intersect. The initial state vector |^q) corresponds to a 
point ipo £ CP™. Therefore the join of ipo and D k is a hyperplane Q k+1 of dimension k + 1 which 
intersects the hyperplane £) Tl_fc_1 at a single point Now take the join of ipo and ipi- The resulting 
line clearly lies in hyperplane Q k+1 and thus hits the hyperplane D k at a single point, and this point 
is the Liiders state fa. The interpretation of ipi, on the other hand, is as follows: if a measurement is 
made to determine simply whether the energy is Ei or not, then in the event of a negative result the 
new state of the system will be the point ipi. 
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Our first requirement will be that {fit} and {a t } must be chosen such that the 
normalisation of \i/) t ) is preserved for all t. We note that the conjugate of ()65|) is 

dty t \ = & t \pldt + (iJ>t\oidW t . (66) 
By virtue of the Ito product rule we have 

d(V t |^> = (d(^IM) + (V»t|(d|Vt)) + (d(^|)(d|Vt)), (67) 
and thus by use of (|63j) and ([SSI) we obtain 

^|t^ = (</4 + At)* + $*t>t) dt + (at + a t ) t dWt, (68) 
where for brevity we use the convenient notation 

{Xt)t - m 

for the expectation value at time t of any operator process Therefore, the 

normalisation condition for \ip t ) is ensured if the operators fl t + fit + a[at and a\ + a t 
have vanishing expectation values with respect to \ipt)- It is a straightforward exercise 
to verify that the most general expressions for the drift fx and the volatility a satisfying 
these conditions are 

fit = -iff, - \a\a t + J t - (Jt)t, (70) 

and 

a t = ik t +L t - (L t )t, (71) 

where {Ht}, {Jt}, {Kt}, and {L t } are arbitrary Hermitian operator-valued processes. 

It should be evident that the process Q is obtained if (a) we let K t and Jt vanish 
for all t, (b) we let L t = ^aHt for all t, where a is a parameter, and (c) we assume 
that {Ht} is time independent. Let us investigate therefore the nature of the additional 
physical conditions that we need to impose on the general norm-preserving dynamics 
in order to ensure that the energy-based model is obtained in accordance with these 
specifications. 

The general norm-preserving model contains the four operator-valued processes 
{H t }, {Jt}, {Kt}, and {L t }. One can think of these operators as representing properties 
of the physical environment in which the quantum system exists. In general, the 
environment is changing in a random time-dependent manner. We shall make the 
simplifying assumption of a 'stationary' environment so that {H t }, {Jt}, {Kt}, and 
{L t } are now replaced by time-independent operators H, J, K, and L. Thus our first 
assumption is that the environment is in a state of stationary equilibrium. 

We give the operator H the usual interpretation as representing the total energy 
of the system. This is justified by the fact that if K, J and L are set to zero then the 
conventional Schrodinger equation is recovered. Next we make an assumption that might 
be called the 'universality of the Hamiltonian'. This is based on the observation that 
the Hamiltonian is the only observable that must exist as an element of the dynamics of 
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a quantum system. If our dynamical law is to be universally applicable to any quantum 
system, then the only observable that can enter the discussion is H, and thus we must 
require that K, J and L are all functions of the Hamiltonian. Thus, in effect, we are 
asking that the system should act as its own environment. It is with this assumption 
that an element of nonlinearity enters the dynamics. 

Our final physical requirement is that energy should be conserved in some suitable 
sense. Now in ordinary quantum mechanics with a time-independent Hamiltonian, the 
expectation value of the Hamiltonian is a constant of the motion. This relation is usually 
interpreted by physicists with a certain looseness of language to mean 'conservation of 
energy', but what it means really is conservation of the expectation value of the energy. 
In a situation where the state is undergoing random changes, the expectation value of 
the energy will also change randomly. We can, nonetheless, impose a slightly weaker 
condition of conservation appropriate to this situation by requiring that the process 
{H t } < t<oo should satisfy the martingale relation 

E[H t \F,] = H, (72) 

for s < t. This relation states that the conditional expectation of the expectation value 
of the energy at time t, with respect to the a-algebra JF S , is the expectation value of 
the energy at time s. This is the sense in which the martingale relation provides a 
characterisation of the principle of energy conservation. 

It follows therefore that we need to analyse the process {H t }, defined as in (jlj), and 
require that its drift should vanish. This will ensure that {Ht} is a martingale, and that 
energy is conserved. In fact, we shall impose a somewhat stronger condition. Let f(x) 
denote a bounded function, and write / = f(H). We shall require that for any such 
operator / the corresponding expectation- value process {ft} defined by 

ft ~ ¥Mh) (?3) 
should be an {^j-martingale. This corresponds to the requirement that not only is the 
energy conserved in the sense discussed above but so is the observable associated with 
any function of the energy. With this condition in place we have a suitably general and 
robust representation of the principle of energy conservation. If we take the stochastic 
differential of f t in (J75J) . then by use of the Ito calculus we find that 

df t = 2 ((J/) t - (J)t(f)t) dt + 2 (<L/) t - (L)t(f) t ) dW t . (74) 

Now the martingale condition on {f t } implies that the drift of {f t } in (J74j) must vanish. 
In other words, we require that the covariance of the two operators J and / should 
vanish for any choice of the function f(x). Thus in particular if we set / = J then it 
follows that the uncertainty of J must vanish in the state \ipt), an d hence without loss 
of generality we may assume that J is a constant multiple of the identity matrix, and 
therefore drops out of the dynamics. 

Finally we consider the roles of K and L in the expression for a t in (|55j). To this 
end we shall examine the dynamics of the squared uncertainty of the operator H in 
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the state \ipt)- Now so far we have through our physical considerations specialised the 
general dynamics ()65|) to the particular case 

d|^> = (-i# - §*fo) \^ t )dt + aMdWt, (75) 

where 

a t = iK + L - (L) t , (76) 

and H, K, and L are time-independent and Hermitian, with the further provision that 
K and L are both given by functions of H. We shall call ()75|) the general stationary 
energy-based stochastic Schrodinger equation. 

Let us therefore investigate the extent to which the general stationary energy-based 
dynamics (fT5j) necessarily leads to state reduction. Writing 

V t = (H 2 ) t - (H) 2 (77) 

for the variance of H with respect to the state \ipt), we obtain 

dV t = d(H 2 ) t - 2{H) t d{H) t - (d(H) t ) 2 . (78) 

Now making use of the fact that both {{H 2 ) t } and {(H) t } are martingales (cf. [5J, we 
immediately infer that {V t } is a supermartingale. In fact, a calculation gives 

dV t = -a((H- (H) t )(L - (L) t )) t 2 dt + 2({H- (H) t ) 2 (L - (L) t )^ dW t . (79) 

It is apparent from ()79|) that the drift of the energy variance process is negative 
in the general stochastic extension of the Schrodinger equation given by (|75j) . This 
demonstrates that the presence of some element of state reduction or relaxation is a 
generic feature of the dynamics of (|75|) . regardless of the specific choice of the functions 
determining K and L. 

This result offers some support to the proposal put forward in Ref. [2j that dynamic 
reduction in quantum theory might be an 'emergent' phenomenon. 

In fact, a short calculation establishes that under the general energy-based dynamics 
(JZSJ), the variance process {V t L } associated with the operator L, defined by V t L = 
({L — (L) t ) 2 )t, satisfies the conditions of being a potential, and admits the following 
representation clS cL conditional variance: 

V t L = EtKLoo - L 4 ) 2 ], (80) 

where L OT denotes the terminal limiting value of the martingale {L t } defined by 
L t = (L) t . 

Thus, provided the eigenstates of L are also eigenstates of H , then (|7Hj) necessarily 
implies a reduction to energy eigenstates. In what follows we shall therefore make the 
simplest choice that ensures this condition, namely, K — and L = \o~H, where a 
is a parameter. Nevertheless, we see that in a general setting there is scope for some 
variation in the dynamics of the state vector from that appearing in ©. In particular, 
we can also consider dropping the stationarity condition. Later in this paper we present 
a useful example of a nonstationary dynamical law. 
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8. Conditional probabilities for reduction 

An important special case of the situation described in ^Hl arises when the observable 
G corresponds to the projection operator IT; onto the subspace of states with energy E { 
(cf. |H IE])- I n this case we have the relations HYii = ELif, ILIL = Sylli, = 1, 

and Eifli = H. The spectrum of H may or may not be degenerate. 

Writing n it = (fli)t for the expectation value of the operator fl, in the state \ip t ) we 
deduce as a consequence of the results of ^Sl that 

d7r it = an tt {Ei - H t )dW t , (81) 

and that 

dv it = -o 2 txI{E { - H t fdt + an it (l - 2n it )(Ei - H t )dW t . (82) 

Here denotes the variance of the operator flj in the state \4>t)- We note that in the 
case of a projection operator the variance takes the simple form 

v it = K it (l —'Ku). (83) 

The random variable k^ has the interpretation of being the conditional probability that 
reduction to a state with energy E{ will occur. In particular, the initial quantities 
K{ = 7r,:o are the Drrac transition probabilities from the initial state l^o) to a state with 
energy 

It is evident that the process {ku} is a martingale, and that this martingale is 
closed by the random variable 

Kioo = 1{H X =E Z }- (84) 

That is to say, 

■K it = E< [K ioo ] . (85) 

As a consequence we see that vu can be written in the form 

v it = Et [{K ioo -E t [n loo ]) 2 } . (86) 

In other words, can be interpreted as the conditional variance of the indicator function 
for collapse to a state of energy Ei. Equation (jHEJ) follows immediately from (jHSJ) if we 
make use of the fact that the terminal indicator function for the energy E{ satisfies 
(^ioo) 2 = K ioo . We see therefore that {v it } is a potential. 

Now we proceed to derive another expression for {7^} that will play a key role in 
the developments that follow. It is well known from the theory of stochastic differential 
equations that an equation of the form (|%Tj) can be integrated. If a positive process 
{X t } satisfies an equation of the form 

dX t = a t X t dW u (87) 

and if f* a 2 ds < 00 almost surely for all t £ [0, 00), then we can write 

X t = X exp( [ a u dW u -\f a 2 u du] , (88) 
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where X is the initial value of the process. If {a t } itself depends in some way on {X t } 
then one cannot say that (|8*8*j) 'solves' (|57J). In that situation (jHHjl should be regarded as 
an integral representation of the stochastic diferential equation (|87jl . Nevertheless we 
may be able to extract useful information about the process {X t } by expressing it in 
this form. In the present case we can integrate (jHTjl to obtain 

Ha = ^ exp (a jf {E { - H u )dW u - \a 2 ^ {E { - H u ) 2 du^J . (89) 

After some straightforward algebraic rearrangement this can be put in the form 



7Tj exp 

H it = 



aE t (W t + a f* H u du) - \a 2 E 2 t 



(90) 



exp (a J* H u dW u + \a* J* H 2 du) 

A further simplification is then achieved if we introduce the {^}-adapted process 
{Ct}o<t<oo defined by the relation 

£t = a [ H u du + W t . (91) 



o 



The process {£t} is evidently a Brownian motion with drift. Making use of the relation 

d£ t = <TH t dt + dWt, (92) 

we can then put ()90)1 in the form 

_ ?Ti exp (aE& - \a 2 E}t) 
Hit — t " " \~- (yd) 

exp [a f Q H s d£ s - \a 2 J Q H 2 ds) 

Finally, we note that since Ei^a — 1> equation (}9"3j) leads us to the following identity: 
exp (o ^ H u d£ u - \ ^ Hldu^J = m ex P {° E & - \o 2 E 2 t) . (94) 
Inserting this relation in ()93|) we obtain the following result. 

Proposition 4. The conditional probability process {hu} for reduction to a state of 
energy Ej takes the form 

7Cj exp (oEjtt - \o- 2 E 2 t) 

where ^t — Wt + cr J * H u du. The energy expectation process {H t } is given by 

Ei^iexp (aE&-\o*E%t) 
1 Z^^P(vE&-^E?t) ' iyDJ 
and the energy variance process {V t } is given by 

Vt = Ei *i(Ei ~ H t ) 2 exp (aE& - \a 2 E 2 i) 
J2i ^ exp {aE£ t - \cr 2 E 2 t) 



We observe, incidentally, that it suffices to specify the value of £ t to determine ir itj 
H t , and V t . In other words, the random behaviour of these quantities is specified entirely 
through their dependence on £ t . 
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9. Information theoretic interpretation of the reduction process 

Given the conditional probability ttu for reduction to an energy eigenstate with energy 
Ei, we can consider the associated information entropy S t . Since the conditional 
probability approaches the indicator function (J%1|) asymptotically, we expect the 
associated entropy to decrease on average. This idea can be put into more precise 
terms as follows: 

Proposition 5. The Shannon entropy process {S t } associated with the conditional 
probability process {itu} is a potential. 

The Shannon entropy (or information entropy) associated with 7t it is defined by the 
expression 

S t = -J^itlniTit. (98) 

i 

This entropy is associated in a natural way with the random density matrix process 
defined by the conditional expectation of the terminal state of the system: 



l^ooX^oo I • (99) 



For each value of t, clearly R t is positive semi-definite and has unit trace. It should also 
be evident that 

S t = -tr{R t hxRt). (100) 
We note that the process {Rt} is distinct from the process {p t } defined by 

j 5 t = E[|Vt><Vt|], (101) 

which is deterministic in t. Thus the state Rt represents the best conditional estimate of 
the terminal state of the system, whereas the state p t represents the initial unconditional 
expectation of the state that the system will be in at time t. Evidently we have Rq = Poo- 
For clarity let us call pt the von Neumann state, and Rt the Shannon state. 

Now if the initial state of the system is a pure state \ipo), with minimum von 
Neumann entropy, then as the reduction proceeds the von Neumann state evolves 
into a mixed state p t with higher entropy. Therefore, the von Neumann entropy 
— tr(p t In pt) associated with the mixed state pt increases from zero to the terminal value 
— tr(poo In poo) = — ^71^^2 7^. On the other hand, the entropy of the initial Shannon 
state Rq is — £\ 7Tj m7Tj, and the entropy of the terminal Shannon state Roo is zero. 

Thus the evolution {p t } of the von Neumann state describes the increase in 
ignorance that results in the statistical description of the system as time moves forward; 
whereas the evolution {Rt} of the Shannon state describes the increase in information 
that results as the measurement outcome is revealed. To put the matter another way, 
the entropy St associated with the Shannon state R t is the negative of the information 
content generated by the information flow {J- s }o< s <t up to time t. In particular, we 
expect the Shannon entropy S t to decrease on average, because more information is 
gained as the collapse process progresses. Finally, when the state has reached an 
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eigenstate, becomes a pure state. The proposition above, which we now proceed to 
prove, asserts that this is indeed the case. 

Proof of Proposition 5. To begin, we determine the dynamical equation satisfied 
by the entropy process {S t }. We note that, as a consequence of the dynamical equation 
(18 lj) for the conditional probability process, and the use of Ito's lemma, we have 

d(\n7r tt ) = -\a 2 {E t - H t ) 2 dt + a(E, - H t )dW t . (102) 

It follows, by another application of Ito's lemma, that 



dS t = -\a 2 V t dt - a ^EiTTulniru - H t s)j dW t , 



(103) 



where {H t } is the energy process and {Vt} is the energy variance process. We observe 
that the volatility of {S t } is the covariance of the energy and the logarithm of the 
conditional probability. Since the drift of {St} is strictly negative we see that the 
entropy process is a supermartingale. 

To show that {St} is a potential we need to show that lim^oo E[SV| vanishes. 
Because the conditional probabilities {ttu} are bounded in the range < < 1, 
the entropy is positive, and is also bounded, and thus lim^oo K[S t ] = E [lim^oo S t ] by 
virtue of the bounded convergence theorem. On the other hand, ()84j) implies that 7Tjoo is 
unity if the terminal energy is Ei and zero otherwise. Therefore, lim^oo S t = almost 
surely, and that establishes the result. {> 

The fact that {S t } is a potential leads to the following observation concerning the 
Shannon entropy and energy fluctuations during the reduction process. 

Proposition 6. The Shannon entropy process {St} is given by the conditional 
expectation of the integrated future energy fluctuation level: 



S t = VE 



2 



VAs 



(104) 



Proof. To derive this result we integrate the dynamical equation ()103|) satisfied by 
the entropy to deduce that 



St — Si 



\a 2 jT V s ds- a jf (j^EiiTislnns-HsSsJ dW s . (105) 

Taking the conditional expectation of this relation we infer, after some rearrangement 
of terms, that 



E t [S T ] = St-±o- 2 E t 

't 



VAs 



(106) 



The identity (|104|) then follows from the fact that {S t } is a potential. {> 
It is interesting to note, incidentally, that if we let t — > in ()106|) . we obtain the 
following formula for the cumulative energy fluctuation during the collapse process: 



¥ 2 



POO 

/ E\V.]ds = S . (107) 
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If the information entropy associated with the initial transition probabilities {7^} is large, 
so that the initial pure state |?/>o) is a highly homogenised superposition of the energy 
eigenstates, one would expect the energy fluctuations during the reduction process to be 
large. Conversely, if this entropy is small, so that the initial state is close to one or a few 
of the eigenstates, then the energy fluctuations during the reduction process should be 
small. Proposition 6 makes this intuition precise. In particular, the right-hand side of 
(JTUZj) measures the entropic uncertainty of the initial energy dispersion (see, e.g., [E]), 
and is independent of the energy spectrum of the system. 

10. Remarks on the ancillary linear dynamics for the state vector 

One of the main goals of this paper is to present in detail a general method for obtaining 
the solution to the dynamical equation (jH)). Before embarking upon this, however, we 
shall first consider the properties of the linear stochastic differential equation 

d|&) = -iH\(f> t )dt - \a 2 H 2 \<p t )dt + §<rff |&>d&, (108) 

and study the relation of this equation to (|HJ). The stochastic differential equation (jl()8|) . 
which we shall call the ancillary equation, plays an important role in the analysis of (j3J) • 
In this section we shall also introduce some change-of-measure formulae that will be 
applied in later sections of the paper. 

In the analysis that follows in this section, the process {|0t)}o<t<oo is to be 
understood as defined on a fixed probability space (f2, .F, Q) with filtration {.Fi}o<t<oo 
with respect to which {£t}o<t<oo is a standard Brownian motion. The precise relation 
of the measure Q appearing here to the measure P introduced earlier will be specified 
shortly, as will the relation between the processes {£ t } and {Wt}. In particular, the 
process {^} introduced in this section has no a priori relation to the process {^} with 
the same name introduced in ^jHl though in what follows the connection between these 
processes is made precise. 

The solution to the ancillary equation (jl()8|) is given by 

|0 i )=e- i ^ + ^-^ 2 ^|0o), (109) 

where |0o) is a prescribed initial state, normalised to unity. The fact that (jlU9j) implies 
(|1(J8|) can be verified by a direct application of Ito's lemma 

d\<k) = mt,t))dt + |0'(&,*)>d& + ||0"(6,t))(d6) 2 , (HO) 
with \(j) t ) = |0(£ t ,;£)), where the function |0(£,£)) is defined by 

|0(e,t))=e-^ + ^-i^ 2t |0 o ). (HI) 

The dot and prime in denote differentiation with respect to t and £, respectively. 

We see that as a consequence of ()109|) that the squared norm of \4> t ) takes the form 

= <0o|e<^-^ 2 U>- (112) 
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Now writing H = £\ E^Yl^ where IT as before denotes the projection operator onto the 
Hilbert subspace for which H takes the value E^ we have 

(<j )t \<j )t ) = J2^ aEA - ha2E?t - (H3) 

% 

In other words, ((j)t\4>t) can be expressed as a weighted sum of geometric Brownian 
motions. Here 7Tj as before signifies the Dirac transition probability 

(114) 

mm) 

from the initial state |0o) to the Liiders state \<pi) associated with the initial state |0o) 
and the eigenvalue Ei. It follows immediately by virtue of the properties of geometric 
Brownian motion that the process {{4>t\4>t)} is a martingale in the Q- measure, satisfying 

E?[<&|&)] = (0 S |0 S ). (115) 

This result can also be seen to follow directly from ()112J) . by an application of Ito's 
lemma, which shows the squared norm of \<pt) satisfies the dynamical equation 

d<&|&) = aH t {<j> t \4> t )^ t , (116) 

where the process {H t } is defined by 

(H<? Ht > 1°™^) 

The stochastic differential equation (jllfij) can then be put in an integral form to give a 
useful alternative expression for the squared norm: 

(&|&) = exp (a H u d£ u - \a 2 jf # u 2 d^ . (118) 

Our intention is to show that the process {H t } defined in (j!17j) can in fact be 
identified with the energy process defined in (J3J). For this, we consider the dynamics of 
the normalised state vector 

\i> t ) = {<f>t\(t>t)- V2 \<t>t). (119) 
If we write N t = (4>t\4>t) 1 ^ 2 for the normalisation factor, then by Ito's lemma we obtain 

dNf 1 = l<j 2 H 2 N- l &t - \aH t N- l <\£ t . (120) 
Hence for the dynamics of the normalised state \if) t ) = Nf 1 ]^) we have 

d|^> = AT^di^) + (diVT 1 )!^) + (diVr 1 )(d|0 t )), (121) 

and thus 

d|^> = -iH\i(>t)dt - y (h 2 + 2HH t - 3# 2 ) \fa)dt + \{H - H t )\^ t )d£ t . (122) 
This expression can be simplified if we introduce a process {Wt} by the relation 

W t = tt-<r f H u du. (123) 



o 
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Then the dynamics for the normalised state vector {|^)} can be written 

d|Vt) = - iH\^ t )dt - \a\H - H t ) 2 \^ t )dt + \a{H - H t )\^ t )dW t , (124) 

which is identical in form to (JHJ), and this leaves us with the problem of the interpretation 
of the process {Wt}. 

Now is by hypothesis a Q-Brownian motion, so evidently {Wt}, as defined in 
(I123|) . is a Q-Brownian motion with drift. We can, however, find a new measure P with 
respect to which {W t } is a P-Brownian motion. The precise statement is as follows. Let 
us fix a finite time T < oo. Then the relevant change-of- measure density Q- martingale 
{$t}o<j<T appropriate for transforming from Q to P over the time horizon t G [0, T] 
is defined by $ t = (<j>t\<f>t)i or equivalently (|118|) . Thus, if A <G Tt denotes any JF T - 
measurable set, and if denotes expectation with respect to the measure Q, then we 
define the probability of the event A with respect to the measure P by the formula 

F(A) = E Q [$ T 1 A ] . (125) 

The theorem of Girsanov ( 20] |23 I2H|) allows us to infer that if {£ t } is a Q-Brownian 
motion, then the process {W t }o<t<T defined by (1123(1 is a P-Brownian motion over the 
given time horizon. 

We note, incidentally, that if {m t } is any Q-martingale, then the process {M t } <t<T 
defined by M t = m 4 /$ t is a P-martingale. In particular, since the process 

(<f> t \H\<f> t ) = (MHfr**-^**^) 

= ^iEie^-^ 2 ^ (126) 

i 

is a Q-martingale (a sum of geometric Brownian motions is a martingale), it follows 
that the energy process {H t } defined by 

is a P-martingale. Therefore, for any finite time horizon [0, T] the dynamics of ((5(1 
can be reproduced by the following procedure. First, we solve the ancillary equation 
(ll(J8() with the required initial condition. Next, the solution thus obtained is used to 
construct the processes {H t } < t <T, {&t}o<t<T, and {|'0*)}o<*<t- Finally, the change- 
of-measure density martingale {$t} is used to change from the 'ancillary' measure Q 
to the 'physical' measure P, which is used to interpret the statistical properties of the 
dynamics of the quantum system. 

With this information at hand, we can now present another useful characterisation 
of the dynamics of the state vector process. We begin with (j3J) and (j3J), and introduce 
the process by use of the relation ([123(1 . The probability space (ft, J 7 , P) and the 
filtration {Ft} are defined, with respect to which {Wt} is a standard Brownian motion. 
We introduce on this probability space the state- vector process {I 1 !/,)} by writing 

\%) = exp f-iHt - - o-Ht)A |Vo>. (128) 
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Then it should be evident that 

hk> = -7==* ( 129 ) 

and hence that is an unnormalised form of the state vector \ipt)- In fact, so is the 
state vector \<j> t ), but \^f t ) and \4>t) have different norms. 

The significance of the process {l^)} is that this process is identical (modulo 
straightforward minor changes in notation) to the nonunitary evolution introduced 
and used by Pearle [2U 123 for the formulation and analysis of collapse models. In 
particular, equation (2.1) of J2H] is identical to our equation (|128|) above. Pearle [25] 
asserts that (|128|) represents "the most transparent formulation of the energy-based 
collapse model". Although ()128j) does indeed represent a formulation of the model, 
it can hardly be regarded as transparent. The problem is that the definition of {£ t } 
involves \^t), and hence ()128j) is, in effect, no more than an integral representation of 
the nonlinear stochastic differential equation (JHJ). To put the matter differently, whereas 
{|^t)} depends on the probability law of the process {£ t } depends on {(^|^)}; 
this is the content of equation (2.2) of J2HJ. Thus, when in what follows we speak of 
obtaining a 'solution' to @, it should be emphasised that we are not merely seeking 
a 'reformulation' such as that represented by (j!28j) . or a change-of-measure induced 
linearisation. 



11. Observation of the energy in the presence of noise 

We now present a general method for obtaining an explicit solution to the stochastic 
differential equation (JHJ). The method we propose ties in very suggestively with the 
theory of nonlinear filtering as developed for example in [2*Tj . 

Let the probability space (fl, J 7 , P) be given, and let {Q t } be a filtration of T 
with respect to which a standard Brownian motion {B t } is specified, together with 
an independent random variable H. We assume that H is ^o _mea surable, and that 
it takes the values {£?i}i=i,2,-,jv with the probabilities {iTi}i=i t 2,- -,N- If we look ahead 
briefly to the results that will eventually follow, the random variable H will have the 
interpretation of representing the terminal value of the energy after state reduction, 
given the Hamiltonian H and the initial state IV'o)- However, for the moment we assign 
no a priori physical significance to H and {B t }, which are introduced as an ansatz for 
obtaining a solution for (JHJ). 

Now suppose we define a random process {£t}o<t<oo according to the scheme 

& = o-Ht + B t , (130) 

where a is a positive constant. Since our units are such that H — 1, the random variable 
H can be thought of as having units of [T] _1 , and hence a, B t , and £ t all have units of 
[T] a . Later we shall identify a with the parameter appearing in the dynamical equation 
(JSJ), but for the moment we leave its value unspecified. 
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The process {£ t } introduced here has no a priori connection with the process having 
the same name introduced in £JH] Nevertheless, as we proceed it will be indicated in what 
sense these processes can be identified with one another. In the probability measure P, 
the process defined by (\VS0\i is a Brownian motion with a random drift rate aH. 
For each value of t one can think of £t as providing noisy information about the random 
variable H. That is to say, given the value of £ t one can try to infer information about 
the value of H. The presence of the independent noise B t interferes with this process. In 
particular, for small values of t, say those for which t <C a 2 , it is typically the case that 
\B t \/ai > 1/a 2 . This follows from the fact that E[\B t \] = v / 2t/ir. Thus if t < a 2 and if 
\H\ <C 1/cr 2 , then knowledge of £t/crt provides little information about the value of H. 
On the other hand, for large values of t we have £t/cr£ « H. We emphasise that at this 
point in our analysis the interpretation of {£t} is irrelevant, since it is being introduced 
as an ansatz for obtaining the solution to (|3*|) . Nevertheless, it will be worthwhile to 
remark as we proceed on various aspects of the nature of the 'information process' {£*}• 

Let {J-f} denote the filtration generated by {^}. We consider the process 
{if 4 }o<t<oo generated by the conditional expectation 



H t = E 



H\J* 



(131) 



Intuitively, conditioning with respect to the a-algebra T\ means conditioning with 
respect to the outcome of the random trajectory {£ s }o< s <i- Clearly, T\ C Gt since 
knowledge of H together with {B s } < s < t implies knowledge of {£ s }o<s<t, although the 
converse is not the case. 

Proposition 7. The conditional expectation E[i/|jF|] represents the best estimate for 
the value of H given the trajectory of the process {C, s }o<s<t from time up to time t. 

Proof. Consider the problem of finding an J-f -measurable random variable Y t 
that minimises the expected value of the squared deviation of H from Y t , given the 
information Tl. Thus we wish to find a choice of Y t that for each u G Q minimises 

J t = E[(H-Y t ) 2 \tf}. (132) 

Since Y t is assumed to be -measurable, we have 

E[(H - Y t ) 2 \^] = E[H 2 \^} - 2Y t E[H\^} + Y 2 . (133) 

Now setting Y t = E[if|jF^] + Z t , where Z t is any jF|-measurable random variable, we 
find that 



Jt = E 



(H-H t ) 2 \J< +Z 2 , (134) 



where H t = E[iJ|jF^]. Therefore, Jt achieves its minimum if and only if Z t = 0. ^> 
The intuition behind this result is as follows. We can think of H as being a hidden 
variable. Its value is hidden by virtue of the noise process {B t }. The best estimate 
available at time t for the value of H is the process {H t } defined by ()131|) . Our goal now 
is to show that {H t } can be identified with the energy expectation process (jU associated 
with the standard energy-based stochastic extension of the Schrodinger equation. 
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12. Optimal estimation of the energy 

We proceed in this section to calculate the conditional expectation (J131)) to establish 
the following useful result. 

Proposition 8. Let H be a random variable taking the value E± with probability 7Tj 
(i = 1,2, ... ,n), and set £ t = oHt + B t for < t < oo, where a is a constant 
and the Brownian motion {B t } is independent of H. Then the conditional expectation 
H t = ~K[H\J-f] is given by 

= E^exp {aE&-\o*E%t) 
1 J2i Ki exp (aE£ t - l<J 2 Eft) ' 1 ' 

Proof. First, we observe that {^} is a Markov process. To establish that is 
Markovian, we need to show that for all T > t the conditional probability distribution 
of £t given the history {£ s }o< s <i is equal to the conditional probability distribution of 
£t given the value £ t of the process at time t alone. In other words, we need to establish 
the following result: 

Lemma 1. Let £t = o~Ht + B t , where H is a random variable taking the values Ei 
(i = 1,2, ...,N) with probability F(H = Ei) = iti, a is a constant, and {B t } is a 
standard ¥ -Brownian motion, independent of H. Then for all T > t and for all x £ R 
we have 

P(fr<4^)=P(6r<a|&). (136) 



Proof of Lemma 1 . It suffices to show that 

P (6 < x\t„ C : -C, U) = P (6 < x\Q (137) 

for any collection of times t, s, si, s 2 , . . . , such that t > s > Si > s 2 > • ■ ■ > s k > 0. 
Now it is a remarkable property of Brownian motion that for any times t, s, s± satisfying 
t > s > Si > one can show that 

B B 

B t and — — are independent. (138) 

s S\ 

More generally, if s > Si > s 2 > s 3 > 0, we find that 

— — and — — — are independent. (139) 

s si s 2 s 3 

In each case the result stated follows after a calculation of the covariance of the indicated 

variables. Next we note that 

£,s ssi B s B Sl 



S Si S Si 

It follows therefore that 



P (6 < X \^3i Csi) Cs 2 ) ■ ■ ■ ) £sj 



(140) 



p 1 


Zt<x 








1 






S 


Sl 


Sl 


S2 


Sk-1 


B s 






B S2 






Bs k \ 


s 


81 


Sl 


1 

S2 


. . . , 


Sk-1 


Sk J ' 



e 



fc 

Sk 



^,ri-^,^-^,...,r^zi_^ . (mi) 



Quantum noise and stochastic reduction 



29 



However, since £ t and £ s are independent of B 3 /s — B Sl /si, B Sl /s\ — B S2 /s 2 , ■ ■ ., 
B Sk _ 1 /sk-i — B Sk /sk, the desired result of Lemma 1 follows. (} 
Continuing with the proof of Proposition 8, we note next that 



E 



H\J* 



E[H\£ t ]. (142) 

That is to say, rather than conditioning on the a-subalgebra T\ generated by {^ s }o<s<< it 
suffices to condition on £ t alone (conditioning with respect to a random variable means 
conditioning with respect to the cx-algebra generated by that random variable). The 
additional information in {£ s }o<s<t does not allow us to improve the estimate of H once 
£t has been given. This follows from the fact that {^} is Markovian and that 

lim - = aH. (143) 

t^oo t 

To calculate E[if|£ t ], we require a version of the Bayes formula applicable when 
we consider the probability of a discrete random variable conditioned on the value of 
a continuous random variable. In this connection we recall for convenience that for 
discrete random variables A and B that take on the values A — A^ {% = 1, 2, • • • , n) and 
B = Bj (j = 1,2, • • ■ , m) with probabilities and rj, respectively, then we have the 
classical Bayes formula 

„ it, t, n F(A = A i )¥(B = B i \A = A i ) , , „ N 

F(A = AAB = BA = — P 144 

or equivalent ly 

, * * , \ Qi^iB = BAA = Ai) f \ 

F(A = AAB = Bj) = =g v , Jl =-; ' - , 145 

3 Eti^(B = B j \A = A,^ K } 

since the marginal probability for the random variable B can be written 

n 

F(B = Bj) = ^2q { F(B = Bj\A = Ai). (146) 

i=l 

Alternatively, instead of conditioning directly with respect to the event B = Bj we 
can condition with respect to the random variable B, and write 

F(A = Ai)F(B\A = A,) 



F(A = Ai\B) - m 

qA?(B\A = A^ 



Zt 1 qMB\A = A i y 

where F(B) is the random variable that takes the value rj = F(B = Bj) when 
B takes value Bj, and P(S|A = Ai) is the random variable that takes the value 
F(B = Bj\A = Ai) when B takes value Bj. Clearly, 

n 

F(B) = J2<l^(B\A = Ai). (148) 
i=i 

For our purpose we need the analogue of (|147|) applicable in the situation where A 

is a discrete random variable and B is a continuous random variable. In that case 

„iT^ F(A = Ai)p(B\A = Ai) , , \ 

F(A = Ai\B) = -± % ±, (149) 
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or equivalently 

f ^^'- ») ' (150) 

since 

n 

p(B) = J2<liP(B\A = A i ). (151) 

8=1 

Here p(x) denotes the density function of the continuous random variable B, so 

F(B <b)= I p(x)dx, (152) 

J — oo 

and p(x\A = Ai) is the conditional density of B given A — A it so 

P(£ < b\A = At)= ( p{x\A = Ai)dx. (153) 

J — oo 

The random variable p(B), resp. p(B\A = Ai), takes the value pip), resp. p{b\A = Ai), 
when B takes the value b. 

Equation ()150j) is the version of the Bayes formula we require in order to determine 
the conditional expectation ()13H) . In particular, since £ t is a continuous random variable, 
we have 

ra/rr F(H = Ei)p(£ t \H = EA 
F{H = EilCt) = -i ^ 

P(ff = E 8 )p(&|ff = £») 
~ £\P(# = EMUH = Ei) 
KiP(£t\H = Ej) 
EiMZtlH = EtY { ] 

Here p(£t\H = Ei) denotes the conditional density function for the random variable £ t 
given that H = Ei. Since {B t } is a standard Brownian motion in the P-measure, the 
conditional probability density for £ t is Gaussian and is given by 

p(£ t \H = E^ = iexp ~ °m 2 ) • (155) 



It follows from the Bayes law ()154|) that the desired conditional probability for the 
random variable H is given by 

, , 7Tj exp (crEiCt — ha 2 Eft) 
F(H = EM t ) = ' V , — 2 . n o \ ' 156 

Therefore, we deduce that 

H t = E[if |&] 

= J2E i nH = E i \£ t ) 

i 

E a inEi exp {aE& - \a 2 E 2 t) 
Y^i TTi exp (aE^ - \o 2 E 2 t) ' 
That concludes the proof of Proposition 7. (} 



Quantum noise and stochastic reduction 



31 



More generally (see, e.g., [S3j), a similar argument establishes that for any bounded 
function x —>■ f(x) we have 

^nJiE^exp^E^-^Eft) 



E 



f(H) 



n 



J2i n exp {aEiit - la 2 Eft) 
13. Existence of the innovation process 



(158) 



We now proceed to establish the following basic result. 

Proposition 9. Let and {H t } be defined as in Proposition 8. Then the process 
{W t } defined by 

W t = i t -a [ H s ds (159) 
Jo 

is an {T\{ -Brownian motion. 

Proof. Starting with the relation £ t = o~Ht + B t we define {W t } as above, with {H t } 
defined as in (|131|) . To show that {W t } is an {.T-^j-Brownian motion it suffices to show 
that {W t } is an {jF^j-martingale and that (dW^) 2 = dt. First, we shall demonstrate 
that {W t } is an {jF^j-martingale. Letting t <T we deduce that 
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(160) 



by Fubini's theorem. Next, we note that 

cT , rt 



E 



ds 



E 



Ft 



ds + y E 

H s ds+ [ H t ds 
o Jt 

[ H s ds + (T-t)H t . 
Jo 



Ft 



ds 



(161) 



Here we have used the fact that {H t } is an {^f }-martingale. Substituting ()161|) in ()160j) 
we obtain 



E 



n 



otE 



H 



n 



+ E 



B 



T 



F' 



a I H s ds. 
'o 

Finally, we observe that by the tower property of conditional expectation we have 
E 



Bt 


f( 


= E 


E [B T \F t B ,H] 


Fl 


= E 


B t 


Fl 



(162) 
(163) 



Inserting this in (|162j) we obtain 



E 


W T 


Fl 


= otH t + E 


B t 


Fl 


-J 














Jo 



HAs 



Quantum noise and stochastic reduction 32 



E 
E 



(atH + B t ) J*\ - at H s ds 
- 1 Jo 



a / H s ds = W t , (164) 



o 



and this establishes that {Wt} is an }-martingale. Next, we observe that since 

dW t = a(H - H t )dt + dB t , (165) 

it follows at once that (dW t ) 2 = dt. Together with the fact that {W t } is an {J-f}- 
martingale we conclude that {Wt} is an {jF^}-Brownian motion. ^> 
We call {Wt} the innovation process associated with the dynamics of the wave 
function. The significance of the fact that {Wt} is an {jF^}-Brownian motion is that 
the process {£ t } as defined in ()13()j) satisfies a diffusion equation of the form 

d£t = o-H t dt + dW t , (166) 

where H t = H(£ t ,i)- As a result, one can prove that T\ = T^ 1 T \ that is to say, the 
information set generated by {Wt} is equivalent to that generated by {&}. It is the 
innovation process {Wt}, and not the noise {B t }, that 'drives' the dynamics of the state 
vector process in (JHJ). 

Now let |-?/>o) be the initial normalised state vector of the quantum system, and let 
Ilj denote for each value of i the projection operator onto the subspace of Hilbert space 
corresponding to the energy eigenvalue Ei, which may be degenerate. As before, we let 

|&) = Trr^n^o) (167) 
denote the Liiders state corresponding to Ei, and we write 

7T it = F(H = E^t) (168) 
for the process defined by (|156|) . 

Theorem 1. The solution of the stochastic differential equation 

d\^ t ) = -iH\A)M - - H t ) 2 \^t)dt + \a{H - H t M)dW t (169) 

with initial condition \ip ) is given by 

i^) = E e " i ^ /2 i^)- ( 17 °) 

i 

Here \<pi) denotes the Liiders state for the eigenvalue Ei, and 

Tr.exp (aEiZt-^Eft) 
lt Ei exp (aEiit - \a 2 Eft) ' 

where £t = aHt + B t . The random variable H takes the value {Ei} with the probabilities 
{fti}, and {B t } is a Brownian motion independent of H . The process {H t } is defined by 
H t = Yli.Ei^it, and the {F^} -Brownian motion {Wt} is given by Wt = £t — o~ J * H u du. 
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Proof. It is a straightforward exercise to verify that (|17U|) satisfies the stochastic 
differential equation (J3J) with the given initial condition. In particular, by applying Ito's 
lemma to (jl71j) and using the relation (jlfifijl . we can verify that {ttu} satisfies 

dn it = a{Ei - H t )n u dW t . (172) 
Then with another application of Ito's lemma we deduce that 

dvr-/ 2 = -\o-\E, - H t ) 2 nl/ 2 dt + ^(E t - H t )^ 2 dW t , (173) 

and with this relation at hand a short calculation shows that ()170j) satisfies (J2J). ^> 

14. Direct verification of the reductive property 

Thus, summing up, the stochastic equation Q can be solved as follows. We let H be 
a random variable taking values {E{\ with the probabilities {7^} defined by (|64j) . or 
equivalent ly by 

(^o|n#o) , 
1 (Ml*) ' { } 

Letting {B t } denote an independent Brownian motion, we define the process {£t}o<t<oo 
by writing £ t = aHt + B t . The solution of (jSJ) is then given by ()170|) or equivalently 

l^v = E, A' 2 exp (-W + \aE t j t - \o*E?t) |fr) 
(^expiaE^-^Eft)) 1 ^ 

where the {J-|}-Brownian motion {Wt} driving {I'l/'t)} in 1)169)1 is given by (|159)l . with 
{H t } defined as in (fTFTjl . 

The fact that (j!57j) defines a reduction process for the energy can be verified directly 
as follows. Suppose, in a particular realisation of the process {Ht}, the random variable 
H takes the value Ej for some specific choice of the index j. That is to say, we condition 
on the event H = Ei. Substituting £ t = aEjt + B t for the corresponding realisation of 
{H t }, we have 

Ej ftiEi exp (o-EiBt - l^E^Ei - 2E j )t) 



H 



Ei exp (aEiB t - \a 2 E t (Ei - E 3 )t) 
E, TTjEj exp (ajEj - Ej)B t - \a\Ej - E,)H) 

Ei vfi exp (a(Ei - E 3 )B t - \cj\E, - E 3 )H) 
KjEj + Ei(^) - E 3 )B t - \a\E-, - Ej) 2 t) 

^ + Ei(^) 7T* exp - - \° 2 ( E i ~ E i) H ) 



(176) 



However, the martingale {Mijt} defined for i ^ j by 

M ijt = exp (a(Ei - E 3 )B t - \a 2 {Ei - Etft) , (177) 
which appears in ()176|) . has the following property: 

lim P (M ijt > 0) = 0. (178) 

t^oo 



Quantum noise and stochastic reduction 



34 



In other words, {M^} converges to zero for large t with probability one. We note, 
incidentally, that a geometric Brownian motion, i.e. a process of the form 

X t = exp (uB t - \v 2 t) , (179) 

has the property that it converges to unity in expectation but to zero in probability. 
That is to say, lim^oo E[XJ = 1 whereas lim t _^ 00 P(X t > 0) = 0. Since 

Ht = TTv^Tm ' (180) 

we see that {H t } converges to the value Ej with probability one. A similar argument 
immediately shows that if H = Ej then for each value of i we have 

lim n it = l{ i= j}, (181) 

t — >oo 

which shows that converges to the Luders state corresponding to the energy 

eigenvalue j with probability one, in accordance with the results noted in 

The advantage of the expression ()157j) is that {H t } and {iV't)} are expressed 
algebraically in terms of the underlying random variable H and the independent 
Brownian motion {B t }. As a consequence, we can directly investigate and verify various 
properties of the reduction process (JHJ) without having to resort to numerical integration. 



15. Identification of independent noise and energy 

In solving the stochastic equation (J3J) we have introduced in t jlll the idea of filtering, 
that is, estimation of the value of the random variable H, given noisy information about 
H, where the noise is induced by an independent random process {B t }. Although the 
method is useful in obtaining an analytical solution to (J3J), the introduction of these 
random variables might appear artificial, because it is not immediately obvious how 
these variables emerge out of the problem specified by (JHJ). Remarkably, however, it 
turns out that we can derive the quantities introduced in fllll from the ingredients 
specified in (jHJ) and (J3|. The aim of this section is to show how this can be achieved. 
We start with the following result. 

Proposition 10. Let {H t } denote the process defined by Q and {^} the process defined 
by ({HI]) ■ The random variables = lim^oo H s and B t — £ t — atH^ are independent 
for all t. Furthermore, the process {B t } is a standard Brownian motion. 

Proof. We begin by establishing the independence of the random variables B t and 
Hoo. To this end we note that it suffices to verify that the relation 

E [e^+^-J = E [e xBt ] E [e yH °°] , (182) 

holds for all x, y. The verification of this property proceeds as follows. Using the tower 
property of conditional expectation (see §2) we have 

E[e* Bt+vH °°] = E [E [ e xBt +y H °° \ 

= E [E [e^+iv-™*) 11 ™ | ;F t w ]] 

= E [e^ E [ e (y-™t) H ™ | Tf\ } , (183) 
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where we have used the {^-"^}-measur ability of the random variable £ 4 in the last step. 
Let us now consider the conditional expectation E ^y- axt ) H °° | appearing inside 
the brackets in ()183j) . By use of the expression for the conditional probability of 
obtained in (jl56j) we deduce that 

To proceed further in determining the outer expectation in (j 183)1 we make use of 
the following subsidiary result. 

Lemma 2. 



7T» exp (itEi(J - \Ela 2 t) = exp (a jf # s d£ s - |<x 2 jf # 2 d^ _ (lg 



5) 



Proof. Let us write $ 4 = $(£t, i) for the left-hand side of ()185|) . where i) is the 
function of two variables defined by 

*) = J] 7Tj exp (e^a - ±£ 2 a 2 t) . (186) 

i 

Then by Ito's lemma we have 

d$ t = <5> t dt + $Jd& + |$"(d&) 2 , (187) 

where the dot and the prime denote differentiation with respect to t and £, respectively. 
Next we observe that (d£ t ) 2 = dt, that $ t + |<I>" = 0, and that = aH t $ t , the last 
relation following from (|157j) . As a consequence we see that satisfies 

d$t = o-ff t $ t d&. (188) 

Finally we note that the integral representation for this stochastic differential equation, 
with initial condition $o = 1, is given by the right-hand side of ([185)1 . ^> 
The key point is that the right-hand side of (|185|) can be used as change-of-measure 
density. Recall that {W t } is a standard Brownian motion in the measure P. Since {H t } 
is bounded and {jF^j-adapted, it follows by Girsanov's theorem that there exists an 
equivalent probability measure Q such that the process {£ f } defined by 

& = W t + a [ H s ds (189) 
Jo 

is a standard Brownian motion in the Q-measure. We let $t denote the change-of- 
measure density in the right-hand side of (J185j) . Then for any {J^ }-measurable random 
variable X t the conditional expectations in these two probability measures are related 
according to the scheme 

KiXt] = ±f?[®tX t ], and E?[X t ]=$Xl^«]- (190) 

Equipped with these results we proceed to determine the conditional expectation 
(|183|) . In particular, if we substitute ()184|) in ()183|) and use the fact that the denominator 
appearing in the expectation is the change-of-measure density $ t , and hence {£ t } is a 
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standard Brownian motion in the Q-measure, we can apply the second identity in (|19U|) 
to deduce that 



E\e xBt+yH °°} = E { 



16 E 



S2 e (y- axt ) E >-\ E 'i a2t e q [ e ( 3S '+ £; » CT )?*] 

i 

7l . e {y-vxt)Ei-±E^ 2 t ^(x+E.aft 

i 

J^ne^Ae^. (191) 



This establishes the relation ()182|) . and hence that random variables B t and ifoo are 
independent. In addition, bonus the result 

E[e xBt ] = e^ 2 ' (192) 

shows that B t is normally distributed with mean zero and variance t. 

To complete the proof that {B t } is a standard Brownian motion we are required, in 
addition to establishing its normality, to verify that the process {B t } has independent 
increments. Alternatively, it suffices to demonstrate that 

E \^B t +y(B T -B t )^ = E j^St] E J e y(B T -B t )] ( 193 ) 

for any nonzero constants x, y. Using the definition for {B t } and the tower property of 
conditional expectation we can write 

E ^ e xBt+y(B T -B t )j _ ]g [ e (x-y)£t+y£T e -(xot+ya(T-t))H x ^ 

Once again from ()156|) we have 



E [ e -(«^cr-t))JT»|^ = Y.i' K i e 



-xatEi-yaiT-^Ei+^TEia-^Ef^T 
^ T Eia-\El^T 



(195) 



for the inner expectation in (|194|) . Substituting (|195|) in (|194|) and noting the fact that 
the denominator in the expectation is the change-of-measure density $t we deduce, 
after some rearrangement of terms, that 

E \^xB t +y(B T -B t )^ _ n . Q -x^tEi-y<j{T-t)Ei-lEl<j 2 T jgQ ^(x-y^t+iV+Ei^T j 



= e ^*+^(T-t). ( 196 ) 

Here we have made use of the Gaussian property 

E Q [ e «e*+6€r] = exp (fE^[(a& + bfr) 2 ]) 

= exp (|(a 2 + 2ab)t + |6 2 T) (197) 

satisfied by the random variables £t and £t in the Q-measure. The result of (1196)1 

establishes ()193|) . and thus we conclude that the process {B t } is normally distributed 

with zero mean and variance t, and has independent increments. Therefore, {B t } is a 

standard Brownian motion. (} 
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16. Finite-time collapse model 

In the foregoing sections we have investigated the properties of energy-based collapse 
models for which state reduction is achieved asymptotically in time. That is to say, 
although for a suitable choice of the parameter a the state reaches the close vicinity of 
one of the energy eigenstates virtually instantaneously, for a strict collapse for which 
the variance process {V t } vanishes identically, we must take the limit t — > oo. There 
are circumstances, however, in which it might be preferable to formulate a model that 
achieves strict collapse in finite time duration. An example for such a model has been 
proposed recently [13] . In what follows we shall apply the methodologies developed 
above to work out the properties of finite-time collapse models. 

The model that we consider here, which gives rise to a finite-time collapse, is given 
by the following stochastic equation: 

d|^> = -iH\^)dt - § [r^j (H - H t ) 2 \^ t )dt + i^~ t (H - H t )\i> t )dW t . (198) 

We deduce immediately from the discussion in £0that the dynamical law ()198j) preserves 
the norm (iptl^t) of the state, and that the associated energy process {H t } is a 
martingale. In particular, a short calculation making use of the Ito calculus shows 
that the energy process satisfies 

dH t = a t V t dW t , (199) 

where we have defined, for convenience, the deterministic function {a t } by 

* = 7^~ v (200) 

and {V t } is the associated variance process. Note that (|198|) can be obtained from (J3J 
by the substitution a — > a t . Thus, ()198|) contains two freely specifiable parameters, 
namely, a and T. The latter will be identified with the time at which the collapse is 
completed. More generally, we may regard the collapse time T as a random variable 
having some density p(T) defined on the positive real line. Then the collapse time T 
itself becomes random; the analysis of this case will be pursued elsewhere. Here we shall 
treat T as a fixed parameter. 

In order to solve ()198j) we shall make an ansatz analogous to the one introduced 
in ([130)1 . Now from the point of view of filtering theory, the collapse of the state in 
the model Q takes place only asymptotically because the 'noise to signal' ratio, whose 
magnitude is of order y/t/t, vanishes only asymptotically as t — > oo. Therefore, in 
order to achieve a finite time collapse we consider the use of a Brownian bridge as the 
source for the noise. A Brownian bridge with duration T can be regarded as a standard 
Brownian motion constrained to take value zero at time t — and also at time t = T. By 
using a Brownian bridge as the source for the noise, the value of the unknown random 
variable H will be revealed in finite time T, since the contribution of noise vanishes at 
that time. Specifically, the magnitude of noise to signal ratio is given by a/(T — t)/tT, 
which vanishes as t — > T. It remains to be shown that the solution of such a filtering 
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problem corresponds to the solution of a finite time collapse model (|198|) . In what 
follows we shall demonstrate that this is the case. 

We thus consider the information process {£t}o<i<T defined in this case by 

& = crtH + B t - l -B T , (201) 

where a is a constant, {B t } is a standard Brownian motion, and if is a discrete random 
variable taking the values {E{\ with probability {ni}. It is evident from definition that 
£o = and that it/o't = H. The process {/3t}o<t<T defined by the combination 

Pt = B t - l -B T (202) 

is a standard Brownian bridge on the interval t G [0, T] satisfying j3 = and (3t = 0. We 
assume that H and {flt} are independent. It should be evident from the definition 
that a Brownian bridge is normally distributed with mean E[$] = and covariance 

1 , „ „ x „ 1 



Cov\p.,Pt] =E 



B s B t - ^(sB t + tB s )B T + j^stB 2 T 

1 " |) (203) 

for s < t. In deriving ()203|) we have made use of the independent increments property 
satisfied by {B t } to deduce that E[B s B t ] = E[B s (B t - B s + B s )} = s. The Brownian 
bridge, on the other hand, does not possess independent increments. 

Our objective now, as before, is to determine the best estimate for the variable H 
given the information concerning the trajectory {£ u }o<u<t of the process from time 
u = to time u = t < T . In particular, the conclusion of Proposition 7 remains valid 
in the present context: that is to say, the estimate that minimises the quadratic error is 
given by the conditional expectation E[if|jF|]. To calculate this conditional expectation 
we shall make use of the following key result: 

Lemma 3. Let £t = atH + (3 t , where H is a random variable taking the values Ei 
(i — 1, 2, . . . , N) with probability F(H = Ei) = -K i} a is a constant, and {f3 t } is a standard 
F-Brownian bridge on the interval t e [0,T], independent of H . Then {£t}o<t<T is a 
Markov process. 

Proof of Lemma 3. To show that is Markovian, we must show that 

ntt<x\J*)=F(Zt<x\Z.) (204) 

for all x G M. and all s, t such that < s < t < T. It will suffice to verify that 

P (6 < x\£ s , £ S1 , U, ■ ■ ■ , U) = P (6 < x\£.) (205) 

for any times t, s, si, s 2 , . . . , Sk such that T > t > s > si > s 2 > ■ ■ ■ > s k > 0. We 
remark that for any times t, s, s\ satisfying t > s > s\ the random variables (3t and 
f3 s /s — (3 S1 / S\ have vanishing covariance, and thus are independent. More generally, for 
s > si > S2 > S3 the random variables (3 S / s— $ 8 Js\ and j3 S2 / S2 — /3 S3 / $3 are independent. 
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We note that £ s /s — £ Sl /si = (3 s /s — fl Sl / s i- It follows that 
P(6<x|£ s ,£ Sl ,£ S2 ,---,£ s J =w(t t <x 



P & < x 



S 


? 

Sl 




) 

S2 


' Sfc-1 


ft* . 


K 


-1 






«2 


' s fc _ 


-1 


Sfc / ' 





Sk 



rs_r* L rsi _ ™ Z!±± - \ . (206) 

S Si Si 



Since £ s and £ t are independent of (3 s /s- /3 Sl /s 1 , f3 s J s 1 - f3 S2 /s 2 , ■ ■ ■, (3 Sk ^Jsk-i -f3 Sk /s h , 
the desired result follows immediately. {> 
Because is a Markov process, the conditional expectation E[i7|^f] simplifies 
to H t = E [-ff|£t] so that we only need to specify the value £ t of the process at time t, 
and not the entire trajectory {£, u }o< u <t- We shall first establish the following result. 

Proposition 11. Let H be a random variable taking the value Ei with probability Hi 
(i — 1, 2, . . . , n), and set £t = atH + B t — (t/T)Bx for < t < T , where a is a constant 
and the Brownian motion {B t } is independent of H . Then the conditional expectation 
Ht = E[/7|jF^] is given by 

E, T ,ex P ( *- Ej r -y^ ) 

Proof. The conditional expectation H t = E [#|£t] can be expressed in terms of the 
conditional probability as follows: 

H t = y £ / E i F(H = E i \£ t ). (208) 

i 

To determine the conditional probability F(H = Ei\£) we note that according to the 
Bayes formula that we can write 

F(H = Ei\£ t ) = (209) 



where 



- o-tE^T 



p&\H = Ek) = y'^yexp [- ^^ ) • (210) 

Expression (|21U|) follows from the fact that conditional on H = Ei the variable £ t in 
(120 lj) is normally distributed with mean atEi and variance 

E [(B t - (t/T)B T ) 2 ] = t(T - t)/T. (211) 

Putting these together, we deduce ()207j) after some rearrangements of terms. ^> 
From the expression ()207|) we can infer directly the property that H t — ► Ek as 
t — > T, provided we set H = E k . Writing for the conditional energy process 
H t (H = E k ), My = Ei — Ej for the difference of energy eigenvalues, and substituting 
& = £* = atE k + (3 t in (J271TJ) . we obtain 

J2i KiEi exp ^ T i t ' j 
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(212) 



H^k \ T-t 

Observe that for each % the numerator in the exponent in (|212j) approaches a strictly 
negative number — ~a 2 u 2 k tT. Hence as t — > T all the exponential terms are rapidly 
suppressed and we are left with the desired outcome: = E k . 



17. Innovation process for finite-time collapse model 

Let us analyse the properties of the process {H t } in (|207|) more closely. By taking the 
stochastic differential of (J2(J7j) we obtain 



dH t = cr t V t 



T — t (Zt-o-TH t yt + d£ t 



(213) 



where V t = 'E[H 2 \J 7 ^] — H 2 is the conditional variance of the random variable H. Clearly, 
there exists a choice of a process {Wt} defined in terms of {Ht} and such that the 
drift term in the dynamical equation ()213|) can be removed, when expressed in terms of 
{W*}. It remains to be shown that such a process is a Brownian motion that derives 
the dynamics of the state (|198|) . We shall proceed by verifying the following result. 

Proposition 12. The process {Wt} defined by 



W t = J {i s - aTH^ds + 6 



(214) 



is an {Jt} -Brownian motion. 

Proof. First we note that the tower property of conditional expectation shows 

E[B t \J*] = E [E[B t \F t B , H]\j*\ = E s [B s ] , (215) 

where we write E s [— ] = E[— It follows from ()20H) that £ t = E t [£ t ], and hence that 
ft is given by 



& = otH + 1 



T 



E t [B ti 



(216) 



We now proceed to establish that {W t } as defined by ()214|) is an {jFfj-martingale. For 
t < u we have: 



E t [W„] = E t [£ u ]+Et 

= E t [iu} + 



[ ^ s -aTH s )ds 



1 1 



-(£ s - aTH s )ds + 



-s J t T-s 



{E t [£ s )-c-TH t )ds.(217) 



Here we have used the fact that {^} and {H t } are {J-'f }-adapted, and that Ef[i7 s ] = H t 
for t < s. Therefore, 



E t [W n ] = E t 
+ 



u 



auH + B U - -B T 



Wt-it 



T 



-E t 



asH + B S - -B T 



ds - oTH 



T-s 



-ds 
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= auH t + W t -^ t + aH t ds - crTH t ds 

Jt T-s J t T-s 

+ E t [B t ] f -J-fl-l) ds 



T J t T-s \ T 
= W t -^ t + atH t + E t [B t ] M -1 

= W t , (218) 

where in the final step we have made use of the relation ()216|) . This establishes the 
martingale property satisfied by {Wt}. On the other hand, (|214J) implies (dW t ) 2 = 
(d&) 2 , whereas flUPJ implies (d&) 2 = dt. It follows that (dW^) 2 = dt, and this 
establishes the assertion that {W t } is an {jF^}-Brownian motion. (} 
We remark that by substituting (J214j) in (J213j) we obtain the dynamics 

dH t = a t V t dW t (219) 

for the process {H t } given in ()207|) . which shows that {H t } is a martingale. On the other 
hand, by taking the stochastic differential of the energy process H t = (ip t \H\ip t ) / (i/jt\i>t) 
using ()198J) we have obtained ()199|) . To show that ()207|) is indeed the energy process 
associated with the dynamics ()198p we must demonstrate that the two processes labelled 
by {Wt} are identical. In particular, we have the following result. 

Proposition 13. The innovation process {W t } defined in ()214|) is the Brownian motion 
that derives the dynamics of the wave function in ()198j) . 

Proof. The stochastic differential equation ()198|) can be given the following integral 
representation 

\ik) = exp (-iHt - \ J* a 2 (H - H s ) 2 ds + \ J* a s (H - H s )dW s ^J \i/> ). (220) 

This can be expressed more concisely as \ip t ) = U t R t \ip Q ), where 

U t = exp ( - iHt) (221) 
is the usual unitary evolution operator, and 

Rt = exp (± jf a s (H - H s )dW s - \ ^ a 2 (H - H s ) 2 ds^) (222) 

is the 'reduction' operator. The square of R t , which we denote by M t , is an operator- 
valued martingale, given by 

M t = exp (J a s {H - H s )dW s - \ J a 2 {H - H s ) 2 ds^) 

exp ( /* Ha s (dW s + a s H s ds) - \ /* a 2 H 2 ds) 

— \ (223) 

exp ( / * H s a s {dW s + cr s H s ds) - \ £ v*H*ds) ' 

Let us now introduce a modified Brownian motion {W t *} according to 

W t * =W t + [ o- s H s ds, (224) 
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so dWt = dW t + o~ t H t dt. While {W t *} is a drifted Brownian motion in the probability 
measure P, we can construct another probability measure Q in which the process {W^*} 
becomes a standard Brownian motion. Then, because H is constant in time, we can 
write M t in the simple form 

M t = i- exp (h a s dW* - \H 2 a 2 s ds^j , (225) 

where 

$ f = exp a s H s dW: - \ J a 2 H 2 ds^j (226) 

is a positive martingale process. 

Recall that (|198|) preserves the norm of l^o)- Therefore, if we assume initially that 
(V'olV'o) = 1) then it follows that (^ \M t \if} ) = 1 for all t. Thus we deduce from (I223|) 



and (12211 that 



$t = <-0o] exp (H^<T,dW t * - \H 2 j\ 2 s dt^ |Vo>- (227) 



As a consequence we can write 



exp (Hf*a s dW t *-lH 2 fio*.dt 

M t = ) ( , (228) 

(^o| exp [HJq a s dW t * - \m jl a*dt) |-0o> 



which has the effect of isolating the dependence of M t on {H t }. In particular, M t depends 
on {H t } entirely through the modified Brownian motion {K^*}. The process {H t } in 
turn is given by H t = (ipt\H\ipt) / {ipt\ipt) , from which it follows that H t = (ipo\HM t \ipo)- 
Therefore, by use of ()228j) we have 



(Vo|# exp [H f a s dW:-\H 2 jy s ds) |^> 

H t = ^ / , (229) 

(^o| exp (H J*a s dW* - ~H 2 fia*ds) |^o) 



which shows that H t can be expressed in terms of {W t *} and t. This is given by 

E< *iBk exp (E, p a s dW; - \E 2 a 2 ds) 

H t = j> r 1 , (230) 

E< * ^p jl a s dW: - \E 2 jl aids) 

where 7Tj denotes the initial probability that the eigenvalue attained is E^. 

Now if the process {H t } obtained in (|2()7j) is the energy process ()23()j) . then from 
the relation J * a 2 ds = a 2 tT/ (T — t), we deduce, by comparison of ()207|) and f)23U|) . that 

it = {T-t) [ -J— dW: (231) 
JoT-s 

must be satisfied. To show that (j231|) is satisfied, we remark first that the stochastic 
differential of (J214|) is given by 

dW t + oT^—H t dt = Tj^-fidt + d&. (232) 
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On the other hand, the differential form of (|224j) is 

1 



dW t + o-T- 



T-t 

Therefore, by comparing ()232|) and 

1 



H t dt = dWt 



t ■ 



we deduce the relation 



t 



Ctdt + dw; 



(233) 



(234) 



This, however, is the differential form of (|231|) . It follows that the process {H t } obtained 
in (|2(J7|) is the energy process (|23U|) associated with the collapse model (|198|) . In 
particular, the process {W t } defined in (J214|) is the Brownian motion that drives the 
dynamics of the state (|198|) . {> 
The above result also shows that the process {£t} defined by fl2()lj) is itself a 
Brownian bridge in the Q-measure. This follows from the integral representation ()23H) 
above, which shows that in the Q-measure, under which {W^*} is a standard Brownian 
motion, {£ t } is a zero-mean Gaussian process with autocovariance given by 
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(235) 



for s <t. Here we have substituted the integral representation (J231)) in the right-side of 
(123 5J) . applied the Wiener- Ito isometry and used the independent increments property 
of Brownian motion. We shall make use of this result to establish Proposition 14 below. 



18. Reverse-construction for finite-time collapse model 

We have demonstrated in the previous section that the closed-form solution to the 
stochastic equation (|198|) can be obtained by use of a nonlinear filtering methodology, 
in which we have introduced a pair of independent random data H and {A}. Conversely, 
starting from the stochastic equation (|198|) . we can derive the existence of such a pair 
of independent data. We let {H t } be the energy process associated with the collapse 
process (j!98|) . and define the process {^} in terms of the energy process {H t } and the 
Brownian motion {Wt} according to 

£t = (T — t) f -J—(dW s + a s H s ds). (236) 
Jo T-s 

Then we have the following: 

Proposition 14. The random variables Ht and (3 t = & — atH^ are independent for all 
t G [0, T\. Furthermore, the process {ft t } is a Brownian bridge. 
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Proof. For the independence of the random variables Ht and (3 t it suffices to verify 
E[e x(3t+yHT ] = E[e xf3t ] E[e yHT ] (237) 
for arbitrary x, y. Using the tower property of conditional expectation we have 



E\e x(3t+yHT ] 



E 



f^E 



^y-atx)H 7 



6 



(238) 



Let us consider the inner expectation E [ e (y-°tx)H T |£j _ Uging 

expressions ()209|) and 

(I210|) for the conditional probability distribution of the terminal energy Ht we deduce 



E 



^(y-atx)H T 



T-t 



(239) 



where the process is defined in fl226|) . Recall now that {$4} is the density process 
for changing the measure from Q to P. As a consequence, we have 
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f^E 
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E^ 



E 



7Tj e (y - CTte) ^ e 



.(240) 



However, the process appearing in (|24U|) is a Brownian bridge under the Q-measure. 
Therefore, the expectation in ()240|) can be computed by elementary methods and we 
deduce, after some rearrangement of terms, that 

E [e^+y^] = eVEl e^" 2 • (241) 



Here we have used the facts that if g is a zero-mean Gaussian random variable with 
variance 7 2 , then E^ 9 ] = e^ 2x2 ; and that the variance of the Q-Brownian bridge {^} 
is t(T — t)/T. This proves the independence of {/3t} and Ht- The result (j241j) also 
establishes that under P the process {fit} is Gaussian, and has mean zero and variance 
t(T — t)/T. To establish {f3 t } is a Brownian bridge, we must show that for s < t 
the covariance of (3 S and (3t is given by s(T — t)/T. Alternatively, we can analyse the 
moment generating function E[e x/3s+ ^*]. We thus proceed as follows. First, using the 
tower property of conditional expectation we have 

= E [e x ^ +y ^E[e- a{xs+yt)HT \Z t )] 
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7T; C 



(242) 



Now we note 



E l 



exp <^ ±E C 



x£ s + {y + 



exp {| (a; 2 S (l - ±) + (y + ff ) 2 t(l - |) + 2x(y + ff -£))}, (243) 
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from which we observe that the dependence on the energy eigenvalues {E{\ in the 
summand of (|242|) drops out. As a consequence, we obtain 

E [e^+^j = exp {f (x 2 s(l - f ) + y 2 t(l - |) + 2zys(l -£))}. ( 244 ) 

It follows that the covariance of /3 S and /3 t for s < i is given by 



° 2 E[e* A+ **] 



dxdy 



x=y=0 



•s(l-|). (215) 



This establishes the assertion that {/3t} is a P-Brownian bridge. {> 
19. Time-change and Brownian bridge 

The asymptotic collapse model and the finite-time collapse model that have been 
investigated in this paper are in fact related by an elementary time-change. In this 
section we shall demonstrate how a finite-time collapse model can be seen to 'emerge' 
from an asymptotic collapse model, and vice verse, by the use of time-change techniques 
applied to Brownian motion. 

We begin by noting the following property of Brownian motion. Suppose f(s) > 
is a continuous monotonic function over s e [0, T] such that 
-t 

f 2 (s)ds^oo (246) 



o 

as t — > T, where < T < oo. Let r(t) be given by the solution of the equation 

,r(t) 

f(s)ds = t, (247) 

o 

or equivalently, 



f2{T{t))= {%) ' (248) 

and let {B t } be a standard Brownian motion. Then the process {X t } defined by 

nr{t) 

X t = / f(s)dB s (249) 



JO 

is a standard Brownian motion. To verify that {X t } is a Brownian motion it suffices to 
show that the covariance of X s and X t for s < t is given by s. This follows on account 
of the fact that since r(t) is deterministic, ()249j) shows that {X t } is Gaussian with mean 
zero. Using Ito's lemma we then find that the covariance is given by 

E[X s X t ] = [ n f(s)ds (250) 



.7 

for s < t. Therefore, from the defining relation ()247j) we conclude that the covariance 
of X s and X t is indeed s. As an example, consider the function 

M = (251) 
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Clearly f(s) is monotonic and satisfies f(s) — > oo as s — > T. Furthermore, we have 

/ 2 (s)d S = -^. (252) 
T — r 

Substitution of (j251j) in ()247j) shows that the relevant time-change in this example is 

r W = rpf- ( 253 ) 

As a consequence, we find that the process {Bt} defined by 

ft+T T 

B t = j Y^7s dBs (254) 

is a standard Brownian motion. From (|25Hjl we find that r(0) = and that r(t) — > T as 
i — > oo. Therefore, the time-change (|253jl has the effect of 'slowing down' the process. 

We can consider, conversely, a time-change that has the effect of 'speeding up' 
the process. For this we require the following variant of the previous result. Suppose 
f(s) > is a continuous monotonic function over s G [0, oo] such that 
-t 

f 2 (s)ds^T (255) 



as t — > oo, where < T < oo. Let r(t) be given by the solution of the equation (|247|) . 
and let {B t } be a standard Brownian motion. Then the process {Y t } defined by 
rr(t) 

Y t = / f(s)dB s (256) 



is a standard Brownian motion. This result can be verified by studying the covariance 
of Y s and Y t . As an example we consider the function 

/W = (257) 

which clearly satisfies the condition ()255j) . It follows that the relevant time-change in 
this example is 

tT 

r(t) = — t , (258) 



or equivalently, 



t = ^j, ( 259 ) 
and that the process {Y t } defined by 

tT 

f—t T 

Yt= / — Tf;dB s (260) 
Jo s + T 

is a standard Brownian motion. We find from (J258|) that r(0) =0 and that r(t) — >• oo 
as t — > T. Therefore, in this example the time-change has the effect of speeding up the 
clock variable r(t). 
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With these results at hand we now proceed to establish the relationship between 
the finite-time and the asymptotic collapse models studied in this paper. To begin we 
recall the definition 

& = atH + & (261) 

for the information process {£ t } in the finite-time collapse model, and the fact that the 
Brownian bridge process {(3t\ admits the integral representation 



s 



dB s . (262) 



Next, we consider the time-change given by (J258|) and define {r] T } by 

Then substituting (J261j) here and using the integral representation (|262|) we obtain 

nr~ t ^ t ^Br + [^- sA B„ (264) 

on account of the integral representation (|262j) . However, recalling the relations (|259|) 
and ()254|) . we deduce that {r] T } can be expressed in the form: 

r] T = aHT + B T . (265) 

This defines a 'standard' filtering problem associated with the asymptotic collapse 
model, if we regard r as the time variable. In particular, the best estimate for H, 
given the observation {J 7 ?}, is determined by the conditional expectation 

H T = E[H\r] T }. (266) 

Clearly we have the relation H T = H t t T y Furthermore, we have 

dH T = aV T dW T , (267) 

where 

W T = 7] T - a [ H s ds (268) 



Jo 

is a standard Brownian motion. From (I263D and (12681) we deduce that 



-L- £ t ) = a H T dr + dW T = aH t —^—— dt + -J— dW t , (269) 
T-t^J (T-t) 2 T-t K J 

where 

Wjt_ = / dW s . (270) 

T 1 Jo 1 ~ s 

Expanding the left side of ()269j) we deduce 



d& + j^r t (^t - o-H t yt = dW t , (271) 

which, if integrated, reduces to the relation ()214j) . In this manner we find that by taking 
the asymptotic collapse model ()265|) and applying the time-change according to 
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which has the effect of 'speeding up' the process, we recover the finite-time collapse 
model (HEED - 

Conversely, given the model 

& = atH + B t (272) 

that solves the asymptotic collapse model, we may consider the time-change given by 
and define 

T T 

Vt = ■=, it = o-tH + B t . (273) 

' t + T s t + T v ' 

However, because of the relations t = tT / (T — r) and ()254|) we deduce that 

' B t = (T-r) [ T T J—dB s , (274) 



t + T " v 'J T-s 

which is the integral representation for the Brownian bridge process with respect 
to the time variable r. As a consequence we obtain 

Vt = otH + /3 t , (275) 

and we thus recover the model that solves the finite-time collapse process. We thus 
obtain the following conclusion: 

Proposition 15. The finite-time collapse model (|198|) can be obtained from the 
asymptotic collapse model Q by means of the time-change defined in (j258|) . The reverse 
transformation is obtained by the time-change defined in (J253|) . 



20. Discussion 



The result of Proposition 9 demonstrates that given the dynamical equation © for the 
quantum state and the associated energy process (JIJ we can deduce the existence of the 
asymptotic random variable and an independent noise process {B t }. The variable 
ifoo carries the interpretation of a hidden variable in the stochastic quantum theory. 
More precisely, because is {jF^j-measurable, its value can only be determined with 
certainty after the collapse has taken place. The 'quantum noise' process {B t } represents 
the 'disinformation' that hides before the completion of the collapse process. 

Whether the energy-based reduction models considered here suffice to describe 
measurements and relaxation phenomena in general in nonrelativistic quantum 
mechanics remains an open issue. There are attempts, for example, to formulate a 
spontaneous collapse of the wave packet in a localised region in space (see, e.g., |H], for 
a recent work in this area). However, localisation of a particle in a small region in space 
typically requires large energy. Indeed, in a generic measurement-theoretic context one 
requires an infinite amount of energy to confine a particle in a finite region [Hj, and 
hence it may be unphysical to speak of a true 'position measurement'. Quantities such 
as the position or momentum of a particle thus represent what an experimentalist can 
estimate from appropriate energy measurements. 
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This observation is consistent with the point of view put forward by Wiener [31] 
that "under the quantum mechanics, it is impossible to obtain any information giving 
the position or momentum of a particle, much less the two together, without a positive 
effect on the energy of the particle examined . . . Thus all coupling is strictly a coupling 
involving energy . . ." The basic idea is that, in order to measure a physical quantity of a 
system the measurement apparatus must interact with the system, and this is achieved 
in the form of interchange of particles (typically photons or phonons). When these 
particles interact with the measurement apparatus they create some form of excitation 
which then allows the device to estimate quantities of interest. In this regard we can 
take the point of view that the energy-based models are of fundamental importance in 
describing random phenomenon involving quantum systems. 
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